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ABSTRACT 


This  report  describes  the  concluding  work  on  a  series  of  studies  concerned  with 
determining,  understanding,  and  trying  to  predict  the  behavior  of  polymethyl¬ 
methacrylate  (PMMA)  and  polycarbonate  (PC)  under  impact.  In  the  present  work  the 
impact  velocity  was  extended  to  2.1,  125,  and  220  m/s  using  a  4.5  mm  steel  ball  as  the 
projectile.  Because  of  dispersion  effects  the  pulse  detected  using  an  instrumented 
Hopkinson  impact  bar  was  markedly  distorted.  A  new  semi-empirical  procedure  for 
correcting  the  pulse  was  developed  that  overcame  difficulties  encountered  in  the 
theoretical  methods.  Analysis  of  the  data  showed  that  the  deduced  Young's  modulus 
of  PMMA  slowly  increases  with  increasing  impact  velocity.  However,  the  derived 
modulus  of  PC  decreases  by  a  factor  of  five  when  the  impact  velocity  was  increased 
from  2.1  to  220  m/s.  This  decrease  in  stiffness  is  unexpected  and  anomalous,  but  is 
beneficial  to  reducing  the  potential  for  impact  damage.  A  theoretically-based  model 
was  formulated  to  “predict”  the  impact  response  of  the  two  materials  from  their  elas¬ 
tic  modulus,  densification  properties,  and  stress-relaxation  laws.  This  modified 
spring/dashpot  model  allows  the  response  of  the  polymers  to  be  continuously  varied 
between  series-like  and  parallel-like.  The  PMMA  was  found  to  require  a  predominant¬ 
ly  series-like  mode  description,  whereas  the  PC  response  was  more  parallel-like.  The 
instrumented  Hopkinson  impact  bar  method  and  analysis  procedures  have  now  been 
developed  to  the  point  where  they  could  be  extended  to  a  much  wider  range  of  poly¬ 
mers  and  impact  conditions,  to  yield  scientifically  interesting  and  technologically  useful 
information  about  polymer  behavior  under  concentrated,  rapidly  applied  loads. 
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SUMMARY 


This  work  concludes  a  series  of  studies  on  the  impact  response  of  two  representa¬ 
tive  polymers:  polymethylmethacrylate  (TMMA)  and  polycarbonate  (PC).  The  goals 
of  this  present  effort  were  (1)  to  measure  the  impact  pulse  produced  when  these  poly¬ 
mers  are  struck  at  velocities  up  to  the  highest  value  feasible  using  a  compressed  gas 
gun,  (2)  to  develop  a  means  of  analyzing  that  pulse  data  so  as  to  reveal  the  response 
of  the  material  to  the  impact,  and  (3)  to  develop  a  model  of  the  expected  response 
derived  from  other  basic  property  information. 

An  instrumented  Hopkinson  impact  bar  was  used  to  detect  the  impact  force  as  a 
function  of  time.  Careful  consideration  and  study  was  given  to  the  selection  of  the 
bar  material  and  geometry.  Poor  acoustic  matching  with  the  specimen  can  result  in 
multiple  reflections  at  the  test  specimen/bar  interface  which  confuse  and  delay  the  sig¬ 
nal  as  it  enters  the  bar.  Also  the  dependence  of  sound  velocity  on  frequency  (disper¬ 
sion)  in  the  bar  material  can  markedly  distort  the  measured  signal.  After  several  tri¬ 
als,  PMMA  was  selected  as  the  best  compromise  bar  material.  Digitally  recorded  sig¬ 
nals  were  obtained  of  the  pulses  produced  by  impacts  of  a  4.S  mm  steel  ball  against  PC 
and  PMMA  targets  at  velocities  of  nominally  2.1,  125,  and  220  m/s.  Reproducibility 
between  replicate  experiments  was  excellent. 

The  pulse  data  required  computational  processing  because  of  the  distortion  due  to 
dispersion.  The  theoretical  techniques  for  doing  this  were  inadequate.  Identification 
of  the  first  arrival  of  the  impact  signal  at  the  sensor,  and  some  basis  for  defining  its 
initial  shape,  are  essential  to  translating  the  pulse  information  into  such  terms  as  the 
reaction  force  versus  penetration  depth.  Several  attacks  were  made  to  solve  the 
mathematically  difficult  problem  of  computing  the  change  in  shape  of  the  leading  edge 
of  a  pulse  as  it  travels  down  the  bar.  A  semi-empirical,  semi-theoretical  solution  was 
found  by  careful  examination  of  the  data.  Using  this  procedure,  the  experimental 
pulses  were  corrected  for  distortion,  and  the  materials  response  could  be  computed. 

The  results  showed  that  PMMA  behaves  much  as  expected  on  the  basis  of  our 
prior  work  at  much  smaller  penetration  velocities.  The  derived  apparent  Young’s 
modulus  increased  perceptibly  with  increasing  impact  velocity  being  3.6,  4.25,  and 
4.4  GPa  at  2.1,  127,  and  219  m/s  impact  velocities,  respectively.  In  contrast  PC 
behaved  anomalously,  having  computed  modulus  values  of  1.4,  0.4,  and  0.3  at  2.1, 
125,  and  221  m/s,  respectively.  The  decrease  in  modulus  with  increasing  velocity  is 
contrary  to  usual  theoretical  expectation  and  reverses  the  trend  observed  in  our  prior 
work  at  lower  penetration  velocities.  Whereas  PMMA  first  developed  radial  cracks  at 
the  highest  impact  velocities,  PC  remained  visibly  undamaged  except  for  the  impact 
crater. 

A  theoretical  visco-elastic  model  was  produced  for  “predicting”  the  impact  pulse 
from  the  Young’s  modulus,  the  densification  volume  ratios,  and  the  stress  relaxation 
parameters.  The  model  also  provides  for  a  continuously  variable  material  behavior 
ranging  from  a  series-like  to  a  parallel-like  coupling  between  the  elastic  and  viscous 
components  of  a  spring/dashpot  representation.  The  calculated  pulse  for  the  PMMA 
material  could  be  made  to  approximate  the  experimental  response  by  considering  this 
material  to  behave  in  a  series-like  manner.  The  calculated  pulse  for  PC  suggested  a 
more  nearly  parallel-like  behavior,  but  could  not  be  made  to  fit  the  experimental 
results  without  varying  the  type  of  behavior  during  the  course  of  the  impact.  This 
model  is  suggestive  of  the  important  parameters  affecting  the  impact  response.  How¬ 
ever,  it  contains  many  approximations  and  requires  critical  testing. 
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1.  INTRODUCTION 

The  response  of  a  structural  polymer  to  a  concentrated  load  such  as  i.  potentially 
damaging  impact  hit  is  both  practically  and  technically  important.  This  report 
represents  the  conclusion  of  a  series  of  studies  of  the  various  phenomena  that  result 
from  impact  encounters  of  polymers  with  small  hard  bodies.  The  long  term  goal  has 
been  to  gain  scientific  understanding  which  in  turn  can  contribute  to:  (1)  the  evolu¬ 
tion  of  polymers  having  improved  impact  resistance,  (2)  the  diagnosis  of  impact- 
related  failure,  and  (3)  improved  techniques  or  criteria  for  characterizing  materials 
with  respect  to  their  impact  resistance. 

The  approach  has  been  to  study  impact  under  conditions  that  closely  resemble  nor¬ 
mal  experience  while  at  the  same  time  trying  to  simplify  and  idealize  the  experiments 
so  as  to  facilitate  fundamentally  meaningful  measurements,  analyses  and  modelling  of 
the  material  processes  occurring  at  the  site  of  the  impact.  Accordingly  the  conceptual 
ideal  has  been  the  impact  of  hard  spheres  against  the  free  surface  of  a  well- 
characterized  semi-infinite  polymeric  solid.  To  approximate  this  ideal  we  have 
focussed  on  two  widely  studied  standard  commercial  structural  polymers,  viz.  polycar¬ 
bonate  (PC),  and  polymethylmethacrylate  (PMMA).  Thick  specimens  of  these  materi¬ 
als  have  been  used  to  simulate  the  semi-infinite  aspects,  and  the  hard  spheres  consist¬ 
ed  of  steel  ball  bearings. 

In  our  preceding  studies(,“6)  these  mateiials  have  been  investigated  at  impact 
velocities  ranging  over  9  1/2  decades  from  0.08  /u,m/s  to  220  m/s.  Covering  such  a 
velocity  span  has  required  several  kinds  of  instrumentation.  Part  of  our  challenge  has 
been  to  establish  concordance  between  the  various  measurement  methods.  Mechani¬ 
cal  test  machines  using  conventional  load  cells  to  measure  the  reaction  force  were  con¬ 
venient  for  penetration  at  constant  velocity.  At  the  other  extreme  free  flight  ballistic 
impacts  were  achieved  using  a  compressed  gas  gun  to  propel  the  projectile.  The  ma¬ 
terial  response  in  this  case  is  sensed  from  the  sonic  impact  pulse  imparted  to  a  Hop- 
kinson  impact  bar.  Translating  the  pulse  information  into  materials  response  informa¬ 
tion  is  much  more  complicated  in  the  case  of  free  flight  impacts  than  for  the  case  of 
slow  encounters  using  mechanical  test  machines.  The  experimental  work  of  the 
present  concluding  study  has  been  devoted  to  the  impact  bar  measurements  over  the 
velocity  range  2.1  to  220  m/s. 

Throughout  this  series  the  observations  of  greatest  interest  have  included  (1)  the 
material  reaction  force  as  a  function  of  penetration  depth  and  time,  (2)  the  stress  re¬ 
laxation  and  strain  recovery  processes  occurring  during  and  upon  halting  the  penetra¬ 
tion,  (3)  the  relative  contributions  to  the  deformation  processes  by  elastic,  anelastic, 
and  viscous  responses,  (4)  densification  beneath  the  indentation  crater,  and  (5)  crack¬ 
ing.  Finally,  a  technical  objective  has  been  to  synthesize  a  model  that  “explains'*  the 
impact  response  in  basic  terms  that  reflect  structural  transformation  and  transport 
processes  in  the  polymers. 

Because  this  report  is  the  final  one  in  this  series,  a  somewhat  extended  review  of 
the  prior  work  is  included  to  provide  perspective  of  the  long  range  technical  goals, 
challenges,  and  concerns  as  well  as  of  the  major  accomplishments,  conclusions,  and  in¬ 
sights.  This  review  also  provides  a  background  for  the  present  studies. 

1.1  Review  of  Prior  Work 

Originally  we  expected  PC  and  PMMA  to  behave  like  linear  visco-elastic  materials 
at  small  strains.  Later,  our  experimental  results  indicated  either  that  linear  visco¬ 
elasticity  model  does  not  apply,  or  that  its  range  of  usefulness  is  imperceptiveiy  nar- 
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row.  Other  results  indicated  a  much  richer  set  of  phenomena  are  involved.  At  the 
start  the  work  was  done  mainly  by  driving  the  spherical  bodies  into  the  polymers  at 
constant  speed  over  the  velocity  range  accessible  through  the  use  of  screw-driven  and 
hydraulic  mechanical  test  machines.  This  was  followed  by  a  completely  different  kind 
of  measurement  technique  in  which  the  spheres  were  propelled  against  the  polymers 
and  the  response  measurements  were  made  using  Hopkinson  impact  b*r  methods. 
This  introduced  an  entirely  new  set  of  problems  that  had  to  be  overcome. 

Ll.l  Prior  Work  Using  Mechanical  Test  Machines 

The  linear  visco-elastic  extension  of  the  classical  Hertz  analysis  for  the  pressing  of 
spherical  elastic  bodies  into  other  elastic  bodies,  provides  a  framework  for  modelling 
the  early  stages  of  the  indentation/impact  process.  The  Hertz  elastic  analysis  gives  the 
following  relationship  for  the  dependence  of  the  penetration  depth  x  on  the  applied 
force  P 


P  -  [4E0Vtf/3(l-  „2)lx^ 


(Eq.  1.1) 


K0x>/i 


where  E0  is  Young's  modulus,  R  is  the  radius  of  the  sphere,  v  is  the  ioisson  ratio, 
and  K0  is  the  quantity  in  brackets.  If  the  material  is  linearly  visco-elastic  it  is  con¬ 
venient  to  introduce  the  concept  of  the  relaxation  modulus  EXt)  defined  as  the  stress 
(at  time  t)  divided  by  a  fixed  strain  applied  instantaneously  at  t  -  0.  For  such  a  ma¬ 
terial  the  equation  corresponding  to  Eq.  1  is: 


Pit)  - 

E0  o  at 


(Eq.  1.2) 


Note  that  if  E{t)  is  time  independent,  then  E(t)  -  so  that  the  two  equations  are 
identical.  It  is  possible  to  invert  Eq.  2  by  the  method  of  Laplace  transforms  as  dis¬ 
cussed  in  Ref.  1,  to  compute  £(f).  At  constant  penetration  velocity  u,  time  can  be  re¬ 
placed  by  x/u.  In  order  to  compute  EXt),  the  force  must  be  measured  as  a  function  of 
x  and  preferable  expressed  analytically  in  terms  that  are  convenient  for  Laplace 
transformation.  The  following  form  was  found(2)  to  describe  the  data  very  well  both  at 
small  and  large  penetration  depths. 


P(x)  -  K0x^ix) 


(Eq.  1.3) 


where 


*(x)  -  11  +  Ax  +  £,e-f>*+  B2(  1  +  ax)<fr>x]/ 
(1  +  B\  +  B2) 


(Eq.  1.4) 


in  which  the  various  constants  are  to  be  determined  from  experimental  data.  The 
justification  for  assuming  linear  visco-elasticity  derived  from  the  fact  that  the  elastic 
Hertz  analysis  shows  that  under  a  ball,  although  the  load  is  concentrated,  the  stress  is 
initially  infinitesimal. 

In  order  to  be  relevant  to  impacts  at  high  velocity,  it  is  necessary  to  determine  EXt) 
for  times  as  short  as  about  1  /is.  This  implies  that  the  load  Pit)  must  reach  a  sensible 
level  in  this  time  interval,  but  this  was  not  possible  using  the  screw-driven  mechanical 
test  machines.  However,  because  the  force-penetration  behavior  was  not  strongly  ve¬ 
locity  dependent,  it  seemed  possible  that  the  short  time  behavior  could  be  deduced  by 
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extrapolation  from  the  data  for  longer  times.  Accordingly,  measurements  were  made 
over  as  wide  a  range  of  indentation  velocities  as  practicable.  An  electro-hydraulic  test 
tnachine(J>3)  extended  the  velocity  maximum  to  about  2.5  m/s  allowing  data  to  be  ob¬ 
tained  over  a  7  1/2  decade  range  of  velocity. 

For  penetration  depths  greater  than  about  1/10  of  the  sphere  radius,  the  observed 
force  increased  essentially  linearly  with  depth  in  all  cases.  Over  this  wide  range  of 
penetration  velocities  there  was  a  150%  increase  in  slope  for  the  PMMA  material  but 
lesis  than  half  of  this  increase  for  PC.  In  fact  for  the  latter  material  the  slope  was  vir¬ 
tually  constant  for  velocities  between  8  nm/s  and  800  nm/s.  Slope  is  related  to 
stillness  and  is  a  measure  of  E{t)  behavior. 

The  premise  of  this  preceding  work  was  that  these  results  could  define  an  effective 
relaxation  modulus.  The  term  “effective”  is  used  because  by  extending  the  computa¬ 
tion  into  the  range  of  large  penetration  depths,  the  visco-elasticity  might  not  remain 
linear.  Nevertheless  this  E{t)  function  would  not  only  serve  as  a  useful  characteriza¬ 
tion  of  the  polymer  but  might  be  used  to  predict  the  impact  response. 

The  variation  of  stiffness  with  increasing  penetration  rate  is  qualitatively  consistent 
with  the  general  behavior  of  relaxation  moduli,  which  are  monotonically  decreasing 
functions  of  time,  i.e.,  the  shorter  the  time,  the  stiffer  the  material  response.  How¬ 
ever,  the  absolute  magnitude  of  the  stiffness  was  several-fold  less  than  that  computed 
from  ultrasonic  measurements  of  the  elastic  modulus  of  PC  and  PMMA.  Hence,  with 
increasing  penetration  velocity  the  stiffness  was  expected  to  increase  to  a  level  con¬ 
sistent  with  the  sonic  results.  However,  unless  the  velocity  dependence  increased 
markedly,  the  measurements  suggested  that  this  would  require  many  more  decades  in¬ 
crease  of  velocity. 

It  gradually  became  clear  that  the  relaxation  modulus  approach  was  inadequate. 
First,  there  was  no  penetration  range  in  which  linear  visco-elastic  behavior  was  experi¬ 
mentally  evident.  Second,  even  with  wide  range  of  experimental  data,  no  clear  limit¬ 
ing  trends  with  further  increase  in  velocity  could  be  detected  in  the  various  parameters 
in  the  <£(x)  function  which  suggest  how  they  would  behave  as  time  approached  zero. 
Third,  the  specimens  after  and  during  the  penetration  process  displayed  a  sharply 
defined  “frozen-in”  densified  zone  that  develops  at  the  site  of  the  impact.  This  zone 
remains  on  unloading.  However,  both  the  zone,  and  the  indentation  crater  disappear 
when  the  polymers  are  heated  to  their  glass-transition  temperatures.  Hence,  the  flow 
that  occurs  during  the  indentation  process  retains  a  memory.  This  behavior  lies  out¬ 
side  of  the  scope  of  a  simple  visco-elastic  model.  Finally,  when  the  rate  of  penetration 
was  abruptly  changed,  unusual  history-dependent  transient  responses  were  observed. 
The  rate  constants  for  the  transients  were  found  to  be  very  nearly  proportional  to  the 
penetration  velocity  just  prior  to  stopping  the  indenter.  For  example,  if  the  penetra¬ 
tion  was  abruptly  stopped,  the  rate  of  load  relaxation  was  found  to  be  proportional  to 
the  velocity  at  which  the  penetration  had  been  produced.  Thus,  the  reason  for  the  ap¬ 
parent  insensitivity  of  the  force-penetration  behavior  to  velocity  during  the  constant 
penetration  velocity  experiments  was  apparently  not  that  the  polymer  was  strain-rate 
insensitive.  Rather,  there  was  an  underlying  marked  strain-rate  sensitivity  which  was 
masked  by  the  simultaneous,  nearly  compensating  stress-relaxation.*3  4) 

The  obseived  oensification*3’  was  about  0.2%  for  the  case  of  PMMA,  and  manifest¬ 
ed  itself  as  an  optically  distinct,  spheroidal  region  having  somewhat  the  appearance  of 
a  glass  bead  embedded  in  the  polymer.  The  m^jor  diameter  of  the  region  was  approxi¬ 
mately  coincident  with  the  bottom  of  the  spherical  indentation.  The  volume  of  the 
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“disturbed"  zone  is  about  6.2  times  the  volume  displaced  by  the  indenter  in  the  case 
of  PC,  and  about  3.6  times  the  volume  for  PMMA.  Interpreting  this  zone  formation 
as  plastic  yielding  lead  to  an  estimated  yield  strain  of  about  0.08  and  0.05  for  PMMA 
and  PC,  respectively,  corresponding  to  yield  stress  of  about  200  MPa  (29  KSl)  for 
both  materials.  However,  other  estimates  lead  to  lower  yield  strains  and  stresses.  By 
following  the  formation  of  the  disturbed  zone  in  PC  using  motion  pictures,  it  was 
found  that  the  zone  develops  almost  from  the  first  instant  of  contact. 

Load  relaxation  as  well  as  creep-recovery  were  studied  at  several  penetration  rates, 
and  penetration  depths  for  both  polymers.  The  load  relaxation  was  found(3)  to  obey  a 
law  having  the  unusual  form: 

Pit)  -  P{0)(\+At)-*  (Eq.  1.5) 

in  which  the  time  t  is  measured  relative  to  the  instant  of  stopping  the  penetration;  A 
and  B  are  experimental  constants,  which  depend  on  the  details  of  the  penetration  his¬ 
tory.  For  both  PMMA  and  PC  the  constant  A  was  found  to  be  nearly  proportional  to 
velocity,  whereas  the  B  term  was  much  more  weakly  dependent  on  this  parameter.  In 
purely  elastic  materials  the  product  AB  must  vanish.  The  data  gave  no  indication  of 
such  a  tendency  with  increasing  penetration  velocity. 

The  above  relaxation  law  is  inconsistent  with  the  concept  of  intrinsic  independent 
molecular  relaxation  processes.  However,  by  formulating  a  relatively  simple  coopera¬ 
tive  molecular  model, it  was  possible  to  account  not  just  for  the  form  of  the  law  but 
also  the  dependence  of  A  on  penetration  velocity.  Briefly  the  model  proposes  that 

(1)  the  fluidity  of  the  polymer  is  increased  by  disentangling  the  polymer  chains, 

(2)  the  disentanglement  process  results  from  the  stretching  or  flow  of  the  polymer  and 
becomes  easier  as  the  unsnarling  progresses,  (3)  the  polymer  spontaneously  reentan¬ 
gles  at  a  rate  proportional  to  its  fluidity.  Thus,  for  any  given  strain  rate  there  should 
be  a  dynamically  equilibrated  state  of  disentanglement,  which  in  turn  Axes  the  fluidity. 
Changing  the  strain  rate  leads  to  transient  responses,  qualitatively  in  accordance  with 
experience. 

In  most  of  the  work  the  standard  ball  diameter  was  4.5  mm.  However,  a  few  mea¬ 
surements  were  made  using  pointed  indenters  of  various  cone  angles.  In  other  experi¬ 
ments  the  use  of  spheres  of  various  sizes  provided  additional  insight  into  the  impact 
process. 

For  a  fixed  penetration  velocity  the  force  can  be  related  to  the  penetration  depth 
and  ball  radius  by  the  empirical  equation: 

P  —  kXmR "  (Eq.  1.6) 


in  which  k  is  a  constant  having  the  dimensions  of  an  elastic  modulus.  Dimensional 
analysis  requires  that  m  +  n  =  2.  As  noted  above  m  is  approximately  1.5  at  small 
depths  but  tends  towards  1  at  larger  depths.  Measurements  at  various  ball  sizes 
showed  n  to  be  0.3-0.4  at  small  depths  and  0.4-0.6  at  large  depths.  At  large  penetra¬ 
tion  depths  the  dimensional  requirements  on  n  and  m  are  violated,  indicating  that  the 
phenomenological  description  is  incomplete.  For  example,  diameter  also  has  a 
significant  effect  on  the  relaxation  processes  through  its  effect  on  the  A  and  B  parame¬ 
ters  in  equation  1.5. 
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1,1.2  Prior  Work  Using  the  Instrumented  Impact  Bar  Methods 

In  the  experiments  at  higher  velocities,  propelling  projectiles  in  free  flight  against 
the  target  materials  adds  an  additional  degree  of  realism  to  the  experimentation.  On 
the  other  hand  the  measurements  become  more  indirect,  since  the  time  dependence 
of  velocity  is  no  longer  independently  controllable.  Load  cells  can  no  longer  be  used 
to  measure  the  reaction  force  when  the  rate  of  load  generates  appreciable  inertial 
forces  in  the  sensor  system.  With  free  impacts  the  pulse  duration  can  typically  be  of 
the  order  of  0.1  ms  or  less.  Under  these  impulse  conditions  measurements  are  gen¬ 
erally  made  instead  using  an  impact  bar.  The  target  material  is  located  on  the  end  of 
the  suspended  bar,  where  it  is  struck  by  a  projectile.  The  bar  is  made  to  be  sufficiently 
long  that  the  resulting  pulse  can  be  completely  generated  and  acoustically  transmitted 
down  the  bar,  before  interference  from  its  reflection  from  the  other  end  of  the  bar. 
The  signal  is  detected  by  a  suitable  strain  gage  or  other  sensor  located  as  close  to  the 
source  as  possible  to  minimize  distortion,  yet  far  enough  distant  from  the  struck  end 
so  that  the  pulse  can  be  assumed  to  be  propagating  as  a  plane  wave.  In  our  case  the 
strain  gages  were  mounted  about  6  bar  diameters  from  the  end.  Measurement  of  the 
impact  momentum,  the  velocity  of  the  projectile  and  if  possible  its  rebound  velocity, 
are  needed  for  subsequent  analysis. 

The  impact  bar  technique  has  long  been  widely  used  to  study  material  behavior  at 
high  strain  rates.  However,  more  conventionally  a  known  shock  or  step  force  is  ini¬ 
tially  applied  uniformly  over  the  entire  cross-section  of  the  specimen.  The  bar  itself 
often  serves  as  the  specimen  and  usually  the  aim  is  to  determine  the  mechanical  prop¬ 
erties  of  the  materials  under  plane  strain  or  plane  stress  condition.  In  the  present 
case,  we  wish  to  study  the  effect  of  a  concentrated  impact  force  of  initially  unknown 
magnitude  or  time  dependence.  Because  the  pulse  is  essentially  generated  by  a  point 
source,  it  propagates  initially  as  a  spherical  wave.  Hence,  in  order  to  analyze  the  pulse 
in  verms  of  plane  waves,  it  is  necessary  to  locate  the  sensors  some  distance  from  the 
specimen. 

If  sound  velocity  were  constant,  independent  of  frequency  or  wavelength,  a  pulse 
would  propagate  without  change  of  shape  and  the  pulse  detected  at  the  sensors  would 
faithfully  reproduce  the  impact  event.  The  pulse  corresponding  to  the  impact  event 
can  be  analyzed  to  compute  the  inferred  penetration  depth,  energy  absorption,  and 
reaction  force.  The  procedure  for  doing  this  is  found  in  the  Appendix.  Thus,  the  im¬ 
pact  bar  data  can  be  converted  to  that  obtained  on  the  mechanical  test  machines. 

However,  sound  velocity  in  a  bar  is  not  constant,  but  depends  on  wavelength  or 
frequency.  The  theory  of  the  velocity  dependence  for  the  case  of  transient  pulses  is 
incomplete.  The  classical  theory  of  dispersion  (velocity  variation)  considers  the  in¬ 
teraction  between  the  dimensions  of  the  bar  radius  and  the  sound  wave  length.  Ma¬ 
terials  having  high  Poisson  ratios  are  more  dispersive  at  high  frequency  than  are  those 
having  low  Poisson  ratios.  However,  in  all  materials  ultimately  the  low  frequency 
components  travel  60%  faster  than  the  high  frequency  ones.  Preliminary  calculations 
show  that  the  range  of  frequency  components  comprising  impact  pulses  unavoidably 
span  the  most  dispersive  region.  It  is  this  kind  of  dispersion  that  is  commonly  dis¬ 
cussed  in  connection  with  impact  bar  methods. 

There  are  other  sources  of  dispersion  at  the  low  frequency  regime.  The  sound  ve¬ 
locity  in  a  bar  is  given  by  y/E0  /p  ,  where  E0  is  the  Young’s  modulus  and  p  is  density. 
However,  in  a  polymer  E0  itself  depends  on  frequency.  (Another  manifestation  of 
this  dependence  is  the  relaxation  modulus).  This  effect  predicts  the  inverse  of  the 
preceding  one,  namely  that  high  frequencies  propagate  faster  than  low  frequency  ones. 
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Both  the  geometric  and  the  molecular  mechanisms  are  usually  considered  in  the  con¬ 
text  of  an  infinite  wave  of  some  fixed  frequency.  These  analyses  forbid  wave  propaga¬ 
tion  in  a  rod  at  velocities  in  excess  of  E0  Ip.  However,  in  pulses  arising  from  a  lo¬ 
calized  encounter,  the  strain  field  is  complex,  and  the  transient  nature  does  not  satisfy 
the  assumption  of  an  infinite  travelling  wave.  The  complicating  transient  effects  are 
discussed  in  2.1,  2.6  and  5.4. 

For  purposes  of  our  studies  if  the  wrong  pulse  shape  is  used  to  compute  the  pene¬ 
tration  vs  depth  behavior,  as  discussed  above,  the  results  can  be  very  misleading.  The 
most  basic  information  about  the  target  materials  is  contained  in  their  response  as 
penetration  first  starts.  Unfortunately,  it  is  the  start  of  the  pulse  that  is  most  severely 
affected  by  dispersion.  In  order  to  ensure  that  corrections  for  these  and  other  factors 
have  been  adequately  considered,  and  because  the  impact  bar  instrumentation  and 
technique  is  less  direct  than  for  the  case  of  mechanical  test  machines,  an  interim  goal 
has  been  to  validate  the  impact  bcr  results.  The  most  straightforward  approach  is  to 
use  an  ideal  elastic  material  such  as  glass  as  the  target  for  which  the  force-time  rela¬ 
tionship  at  the  impact  site  can  be  computed  from  elastic  theory. 

The  two  studies(5,6)  preceding  the  current  one  were  largely  devoted  to  (1)  finding 
suitable  experimental,  theoretical,  and  analytical  approaches  to  overcoming  the  disper¬ 
sion  problem,  or  in  rendering  it  insignificant,  and  (2)  verifying  the  impact  bar  results. 

An  additional  factor  of  concern  is  the  acoustic  coupling  between  the  specimens  and 
the  impact  bar.  Mismatch  of  acoustic  impedance*  results  in  only  partial  transmission 
of  the  signal  from  the  specimen  into  the  bar,  the  remainder  being  reflected  back  into 
the  specimen.  Thus,  large  impedance  differences  trap  the  pulse.  This  also  can  lead  to 
gross  discrepancy  between  the  transmitted  signal  and  the  instantaneous  force-time  re¬ 
lationship  provided  by  the  impact.  The  original  approach  was  to  make  the  impact  bar 
out  of  PMMA,  thereby,  essentially  eliminating  the  acoustic  mismatch  problem.  The 
plan  for  correcting  for  dispersion  was  to  Fourier-analyze  the  observed  pulse  signal 
The  originating  pulse  could  then  be  derived  by  resynthesis  after  correcting  each  term 
for  the  phase  shift  induced  by  the  frequency-dependent  velocity.  The  dependency  of 
the  velocity  on  frequency  would  need  to  be  measured.  However,  the  use  of  a  Fourier 
series  description  was  shown  to  lead  to  artifacts  because  the  pulse  is  non-periodic. 
Furthermore,  correctly  assigning  the  initial  phases  and  carrying  out  the  phase  correc¬ 
tions  becomes  increasingly  uncertain,  difficult,  and  impractical  for  all  except  the  first 
ten  or  so  terms. 

Other  approaches  were  to  use  other  less-dispersive,  elastic  bar  materials.  One  of 
the  additional  problems  with  the  PMMA  bar  was  that  even  though  impact  signal 
strains  in  the  bar  were  very  small,  typically  .0003,  nevertheless,  the  strain  did  not  re¬ 
turn  to  its  initial  unstressed  level  after  passage  of  the  impact  pulse.  This  created 
difficulties  with  respect  to  applying  conservation  of  momentum  calibration  procedures 
to  the  data.  Hence,  the  PMMA  bar  was  replaced  by  an  aluminum  one.  This  had  the 
additional  advantage  that  the  impedance  of  A1  closely  matches  that  for  glass.  The  ob¬ 
served  pulse  obtained  with  the  glass  target  was  very  close  to  the  theoretically  calculat¬ 
ed  pulse.  Results  using  SiC  as  the  target  gave  similar  agreement.  Hence,  the  A1  bar 
was  considered  to  be  performing  satisfactorily,  and  was  used  to  determine  the  impact 
response  of  PC,  PMMA,  and  cured  epoxy  specimens  by  4.5  mm  and  19  mm  diameter 
strikers  at  impact  velocities  of  1.8  and  2.9  m/s. 

*  Acoustic  impedance  Z  is  defined  as  V  E0p ,  and  the  reflection  at  the  boundary  between  materials  1  and 

2  is  given  by  (Z,  -  Z2)/(Z,  +  Z2). 


These  data  were  examined  first  using  such  general  shape  parameters  as  the  peak 
strain  and  the  pulse  width  at  half  the  peak  strain  value.  The  amount  of  energy  ir¬ 
reversibly  consumed  in  the  impact  process  was  found  to  increase  in  the  order  glass 
<  PMMA  <  epoxy  <  PC,  and  that  impacts  involving  large  diameter  projectiles  are 
more  elastic  than  for  the  case  of  small  diameter  ones.  The  data  were  further  analyzed 
to  yield  the  time  vs  penetration  depth  behavior  as  well  as  the  energy  absorbed  in  the 
immediate  region  of  the  impact.  The  force  vs  penetration  depth  results  were  reprodu¬ 
cible,  and  were  similar  to  those  obtained  using  the  electro-hydraulic  test  machine. 
However,  the  absolute  value  of  the  forces  was  20-25%  lower  than  the  value  obtained 
using  the  latter  equipment.  The  reason  for  this  discrepancy  is  not  known. 

The  observed  dispersion  effects  were  smaller  in  the  A1  bar  than  in  the  PMMA  bar. 
Nevertheless  they  became  increasingly  detectable  as  the  impacts  became  more  abrupt. 
In  the  case  of  glass  and  SiC  targets  the  observed  pulse  seemed  to  “taper  on”  over  a 
time  span  of  2.3  /u,s  and  6  /u.s  respectively.  Inasmuch  as  impact  velocities  of 
—  250  m/s  were  planned,  during  this  “taper-on”  interval  the  projectile  would  travel 
about  0.5  -  1.0  mm.  If  this  taper-on  were  interpreted  as  the  actual  force  behavior  it 
would  result  in  a  false  computed  force-penetration  relationship  in  which  the  glass  or 
SiC  was  apparently  very  soft  over  the  depth  range  corresponding  to  this  interval  as  just 
discussed.  This  artifact  will  also  occur  when  the  target  is  polymeric  and  is  impacted  at 
high  velocity.  Therefore,  the  problem  of  dispersion,  while  rendered  innocuous  at  the 
relatively  low  impact  velocities  remained  unresolved  with  respect  to  impact  velocities 
achievable  using  a  compressed  gas  gun.* 

Finally,  in  a  parallel  series  of  experiments  photographs  of  the  impact  of  a  4.5  mm 
steel  projectile  with  PC  and  PMMA  targets  were  made.  Velocities  of  about  100  and 
140  m/s  were  produced  using  the  gas  gun  and  the  pictures  were  taken  with  an 
Ima-Con  high  speed  camera  at  intervals  of  5  /us  per  frame,  and  the  polymer  samples 
were  attached  to  a  massive  steel  block.  In  the  PC  material  a  diffuse  densified  zone 
found  under  the  projectile  during  penetration  and  became  much  more  distinct  after 
the  impact.  No  cracking  occurred.  In  the  case  of  PMMA,  the  zone  was  much  more 
distinct  initially  but  seemed  to  contract  very  substantially  after  impact.  Cracking  in 
this  case  occurred  at  140  m/s  and  initiated  at  the  point  of  maximum  penetration. 

1.2  Objectives  of  the  Present  Study 

The  purposes  of  the  present  contractual  effort  were  (1)  to  conduct  impact  measure¬ 
ments  on  the  PC  and  PMMA  materials  at  the  highest  velocities  feasible  with  the  gas 
gun,  (2)  to  find  a  way  to  correct  the  pulse  signal  so  as  to  reveal  the  actual  force-time 
history  of  the  impact  either  by  rendering  pulse  distortion  negligible,  or  by  making  the 
necessary  corrections,  and  (3)  to  bring  together  the  body  of  information  gained  to  date 
into  a  model  of  the  materials  properties  retevent  to  impact  resistance  per  se  and  that 
distinguished  the  PC  from  PMMA. 

The  prior  work  provided  the  foundation  for  these  continuing  studies.  However,  a 
clear  approach  to  dealing  with  objective  (2)  was  not  evident  although  some  additional 
clues  were  available.  Extending  the  work  into  the  higher  velocity  regime  was  techni¬ 
cally  motivated  by  the  fact  that  in  this  regime  PC  exhibits  continuing  resistance  to  im- 

*  The  problem  of  pulse  shape  distortion  due  to  dispersion  has  long  been  recognized.  In  principle  the 
solution  to  the  problem  is  contained  in  the  classical  theory  of  (time  dependent)  continuum  mechanics. 
Formal  exact  solutions  relating  to  the  present  problem  have  been  given(8'.  These  require  the  evalua¬ 
tion  of  intractable  integrals.  Approximate  solutions  are  known  which  are  only  valid  for  tong  times 
after  the  arrival  of  a  pulse.  These  account  for  the  “ringing”  that  occurs  in  the  trailing  regions.  How¬ 
ever  this  writer  knows  of  no  solution  to  the  behavior  of  leading  edge  as  the  pulse  propagates. 
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pact  fracture,  whereas  PMMA  is  finally  starting  to  exhibit  brittle  cracking.  Hence,  the 
most  pertinent  information  is  likely  to  be  obtained  at  these  velocities.  Furthermore, 
this  range  of  velocity  is  of  most  practical  interest  to  common  experience. 

The  following  section  provides  an  in-depth  overview  of  the  attack  on  these  objec¬ 
tives.  This  is  then  followed  by  detailed  documentation  and  discussion  of  those  ele¬ 
ments  directly  leading  to  the  resolution  of  those  objectives. 
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2.  OVERVIEW  OF  PRESENT  WORK 

Although  our  prior  work  served  as  the  foundation,  the  present  study  has  been 
more  exploratory  than  expected  because  the  means  of  analyzing  the  experimental 
results  had  not  been  resolved.  A  particularly  nagging  problem  that  remained  was  how 
to  identify  the  arrival  time  of  the  leading  edge  of  the  impact  signal,  and  how  to  deter¬ 
mine  the  shape  of  the  start  of  the  pulse.  The  experimental  and  analytical  aspects  of 
the  work  were  interdependent  because  the  choice  of  the  bar  material  and  geometry 
affected  the  strength  and  distortion  of  the  pulse  signal,  which  then  required  correction 
by  the  data  analysis  procedure.  The  path  for  resolving  these  issues  contained  dead 
ends  along  the  way.  This  chapter  provides  an  overview  of  this  study  as  it  proceded, 
including  those  procedures  that  worked,  as  well  as  those  having  a  more  will-o'-the- 
wisp  character.  The  “successful”  portions  are  treated  in  depth  later. 

2.1  Correction  of  Pulse  Shapes  Using  Transfer  Functions 

Our  prior  work  showed  that  a  Fourier  series  approach  to  correcting  a  pulse  slope  for 
dispersic  i  is  fundamentally  inappropriate  and  inadequate.  The  correct  procedure 
would  be  to  use  a  Fourier  integral  representation,  but  this  is  intractable.  However, 
another  approach  widely  used  in  the  signal  communications  field  showed  promise. 
This  method  analyzes  pulses  in  the  amplitude-time  domain  rather  than  in  the 
amplitude-frequency  domain. 

In  the  language  of  signal  communications  an  input  signal  given  by  some  time 
dependence  HO  is  introduced  into  a  “channel”  (amplifier,  transmission  line,  filter,  or 
other  system)  from  which  an  output  signal  G(t)  emerges.  The  properties  of  the  chan¬ 
nel  are  assumed  to  be  fixed  and  not  altered  by  the  signal.  The  channel  can  be  charac¬ 
terized  by  its  own  function  MO,  called  the  transfer  function.  This  concept  can  be 
represented  schematically  as  follows: 

Input  Signal  — 1 *  [Channel]  — *  Output  Signal 

HO  MO  G(t) 

In  our  case  the  input  signal  is  the  impact  event,  the  output  signal  is  that  observed  by 
the  strain  gage  sensor,  and  the  channel  is  the  impact  bar.  The  way  the  channel  acts  on 
the  input  signal  to  yield  the  output  signal  is  given  by: 

Git)  -  f  FU-t)H(t)<Jt  (Eq.  2.1) 

This  equation  can  be  solved  for  Hit)  using  Laplace  transform®  if  Git)  and  Fit)  are 
known.  Since  the  impact  bar  between  the  impact  site  and  the  sensor  location  remains 
constant  for  any  particular  impact  bar,  once  Hit)  for  the  bar  is  determined,  it  should 
be  applicable  to  any  pulse.  Similarly,  by  inverting  Eq.  2.1,  one  can  write: 

Fit)  -  f  G{t-r)H0ir)dr  (Eq.  2.2) 

t'o 

where  H0it)  is  the  transfer  function  for  the  inverse  process  of  correcting  the  output 
signal  to  reveal  the  original  signal.  This  approach  is  widely  used  for  applications  rang¬ 
ing  from  designing  correction  filters  for  high  fidelity  sound  systems,  to  deciphering 
electronically  scrambled  information.  The  transfer  function  approach  is  actually 
mathematically  equivalent  and  related  to  the  Fourier  integral  method.  When  the  out¬ 
put  data  is  produced  at  discrete  intervals,  the  integrals  in  Eqs.  2.1  and  2.2  are  replaced 
by  corresponding  summations. 


For  the  present  case  Fit)  can  be  calculated  on  the  basis  of  the  Hertz  elastic  theory 
provided  that  the  target  material  is  a  purely  elastic  material  such  as  glass  or  silicon  car- 
bide.  The  output  function  Git)  is  the  measured  pulse  as  detected  by  the  sensor  that 
corresponds  to  the  calculated  pulse.  It  is  useful  to  think  of  the  transfer  function  as  a 
kind  of  memory  which  reads  an  incoming  signal  at  some  instant  but  retransmits  it  out 
over  a  period  of  time  with  some  fixed  weighting  pattern.  This  period  is  expected  to  be 
finite  in  duration  ranging  over  the  interval  defined  by  the  maximum  and  minimum 
sound  velocities. 

When  the  start  of  the  observed  pulse  was  used  to  calculate  the  transfer  functions, 
the  results  were  found  to  be  very  sensitive  to  the  value  assumed  for  the  arrival  time 
of  the  leading  edge  of  the  pulse.  Other  determinations  of  H0it)  can  be  obtained  by 
using  other  portions  of  the  pulse.  The  computed  transfer  function  depended  on  the 
portion  of  the  pulse  selected  as  well  as  on  the  value  of  the  period  and  the  leading  edge 
arrival  time.  Noise  and  systematic  variations  between  nominally  identical  experimental 
runs  also  affected  H0it).  To  find  the  best  overall  representation  a  computer  program 
was  developed  that  produced  an  average  H0it)  by  scanning  a  range  of  data.  The 
period  and  pulse  arrival  time  could  also  be  systematically  varied.  That  combination  of 
parameters  that  led  to  a  transfer  function  that  resulted  in  the  least  error  between  the 
calculated  Fit)  and  the  actual  Fit)  was  assumed  to  lead  to  the  best  estimate  of  //„(/). 
The  effort  then  turned  to  obtaining  precise,  high  reliability  data  for  Git ),  for  cases 
where  Fit)  could  be  calculated  with  a  similar  level  of  assurance.  However,  the  accu¬ 
rate  determination  of  the  leading  edge  arrival  time  remained  a  matter  of  high  priority. 

From  the  start  of  the  effort  in  the  transfer  function  approach,  it  was  clear  that  a 
difficulty  would  arise  if  the  dependence  of  Git)  as  a  power  of  /were  to  be  greater  than 
that  of  Fit).  This  emphasized  the  importance  of  the  time  /  which  is  measured  from 
the  start  of  the  input  signal,  i.s.  from  the  arrival  time  of  the  leading  edge.  If  such  a 
power  dependence  were  to  occur  there  could  be  no  mathematical  solution  to  the  con¬ 
volution  equation  2.7.  Therefore,  a  physical  argument  indicates  that  Git)  cannot  have 
such  a  greater  power  dependence  in  as  much  as  Fit)  is  in  fact  the  predecessor  of  Git). 
This  consideration  placed  a  restriction  on  the  acceptable  functional  description  for  the 
leading  edge  portion  of  the  pulse.  However,  this  restriction  was  considered  to  be  help¬ 
ful  rather  than  a  problem  because  it  somewhat  constrained  the  range  of  possible 
values  for  the  start  of  the  pulse  signal. 

As  will  be  discussed  later,  this  resolution  of  the  potential  problem  appears  to  be  in 
error.  Thus,  the  use  of  a  transfer  function  to  correct  the  pulse  for  dispersion  was  later 
abandoned  in  spite  of  its  apparent  theoretical  soundness. 

2.2  Theoretical  Studies  of  the  Behavior  of  the  Pulse  Leading  Edge 

The  need  for  correctly  and  accurately  determining  the  start  of  the  impact  pulse,  be¬ 
cause  spurious  effects  may  result,  has  already  been  emphasized.  The  problem  arises 
because  of  the  difficulty  in  defining  when  and  how  the  pulse,  (which  starts  out  with 
zero  amplitude  and  slope)  rises  from  the  background  noise  level.  As  suggested  by 
Eq.  1.1,  the  original  impact  pulse  is  expected  to  increase  in  force  and  level  as  about  the 
3/2  power  of  impact  time  in  accordance  with  the  Hertz  analysis.  Whether  this  power 
dependence  is  preserved  as  the  pulse  propagates  down  the  bar  was  an  unanswered 
question.  Determination  of  the  start  of  the  transmitted  signal  would  be  greatly 
simplified  if  the  assymptotic  behavior  of  the  leading  edge  could  be  theoretically  estab¬ 
lished.  Accordingly  a  study  was  undertaken  to  determine  the  evolution  of  the  leading 
edge  of  a  pulse  that  initially  is  proportional  to  t3/2.  The  approach  was  to  introduce  this 
initial  pulse  form  into  the  elasticity  equations  and  to  concentrate  on  the  behavior  just 
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at  the  leading  edge.  The  simplest  appropriate  equation  that  incorporates  dispersion 
was  given  by  Love(7)  and  is: 

u„  -  K2v2uxxl!  -  C2^  (Eq.  2.3) 

where  u  is  displacement,  C  is  the  extensional  velocity,  K  is  the  radius  of  gyration,  u 
the  Poisson  ratio,  x  axial  distance  and  rtime. 

Attempts  to  use  this  equation  to  solve  this  problem  analytically  by  power  series, 
Laplace  transform,  or  Fourier  analysis  methods  led  to  either  intractable  or  discourag- 
ingly  complex  formulations  and  were  abandoned.  A  finite  difference  approximation 
was  then  used  to  explore  the  behavior,  assuming  initially  that  the  pulse  propagated 
without  dispersion,  and  then  correcting  the  form  by  successive  approximations  using 
the  above  equation.  As  the  propagation  time  was  made  to  advance,  the  computer 
results  at  the  leading  edge  became  unstable  with  respect  to  successive  iterations.  This 
instability  occurred  independently  of  the  refinements  used  to  represent  the  derivatives 
or  the  scale  of  the  difference  intervals  selected.  It  has  been  already  be  noted  that  for¬ 
mal  solutions*87  to  the  present  problem,  while  statable,  are  veiy  difficult  to  apply  to  the 
leading  edge  behavior. 

The  approach  that  proved  to  be  the  most  fruitful  was  a  semi-empirical  one  based 
on  examination  of  the  experimental  results  in  combination  with  satisfying  certain  basic 
physical  principles.  This  is  discussed  in  Section  5.4  after  discussing  the  analysis  of  the 
experimental  results  in  Section  5.3. 

2.3  Impact  Bar  Dispersion  and  Acoustic  Impedance  Matching  Considerations 

Although  the  previous  study  showed  the  2.5  cm  O.D.  aluminum  impact  bar  to  be 
appreciably  less  dispersive  than  the  PMMA  bar  first  used,  the  A1  bar  was  in  turn  re¬ 
placed  by  a  1.25  cm  O.D.  vitreous  silica  bar.  The  smaller  diameter  was  chosen  to  im¬ 
prove  the  frequency  dependence  of  the  dispersion  and  to  provide  greater  strain  signal 
output.  Vitreous  silica  was  selected  to  reduce  further  the  dispersion  level  because  of 
its  relatively  low  Poisson  ratio.  Furthermore,  Si02  and  Pyrex  glass  specimen  targets 
are  very  well  matched  acoustically.  All  of  these  factors  are  favorable  for  the  determi¬ 
nation  of  the  characteristic  transfer  function  of  the  Si02  bar  material.  However,  the 
acoustic  impedance  match  with  polymeric  target  specimens  is  poor. 

The  choice  of  the  ideal  impact  bar  for  use  with  polymer  target  materials  presents  a 
mild  dilemma.  If  a  polymer  bar  is  used,  the  acoustic  coupling  is  good,  but  the  disper¬ 
sion  effects  are  relatively  severe.  If  a  Si02  or  A1  bar  is  used,  the  acoustic  coupling  is 
poor  (an  increasing  problem  as  the  impact  velocity  increases)  but  the  dispersion  is  less. 
Correction  for  acoustic  mismatch  is  done  by  the  use  of  a  convolution  procedure  as  in 
the  case  of  the  transfer  function  correction.  Both  the  acoustic  impedance  and  the 
acoustic  reverberation  time  as  determined  by  the  sound  velocity  depend  on  the  ap¬ 
propriate  elastic  modulus  and  density  of  the  material. 

The  resolution  to  this  Hobson’s  choice  was  effected  through  the  following  stra- 
tegem.  First,  a  glass  target  mounted  on  the  Si02  rod  was  impacted  to  enable  the 
transfer  function  for  this  rod  to  be  determined.  A  hemisphere  of  glass  was  then  sub¬ 
stituted  for  the  glass  specimen.  The  Si02  impact  bar,  suspended  by  fine  wires  in  the 
manner  of  a  pendulum,  was  mounted  in  tandem  with  a  similarly  suspended  PMMA 
impact  bar,  onto  which  a  PMMA  specimen  was  mounted.  The  Si02  bar  was  then  dis¬ 
placed  axially  away  from  the  other  bar  and  released  so  that  when  it  swung,  the  two 
bars  impacted  end  on.  The  impact  signal  was  read  from  each  bar.  Since  one  impact 
event  generated  the  pulse  in  the  two  bars,  the  two  observed  signals,  when  property 
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corrected,  must  be  identical.  The  pulse  in  the  Si02  bar  was  corrected  by  means  of  its 
previously  determined  transfer  function.  This  corrected  signal,  in  combination  with 
the  observed  signal  in  the  PMMA  bar  now  provided  the  required  information  to  com¬ 
pute  transfer  function  in  the  PMMA  bar.  Thus,  the  complications  due  to  impedance 
mismatch  effects  have  been  eliminated  in  both  bars.  Having  now  (apparently)  achieved 
a  clean  determination  of  the  transfer  function  in  the  PMMA  bar,  that  now  became  the 
preferred  bar.  This  choice  was  also  predicated  on  the  difficulties  that  would  be  experi¬ 
enced  in  correcting  the  data  for  reflection  if  the  Si02  bar  were  used.  The  reverbera¬ 
tion  time  would  he  about  12  /Lts  in  a  12  mm  thick  specimen.  Neither  this  time  nor  the 
reflectivity  would  remain  constant  because  of  the  intrinsic  dispersion  of  the  polymer 
materials.  Since  pulses  of  the  order  of  only  100  ^s  duration  were  anticipated  in  the 
high  velocity  impacts,  these  characteristics  would  be  very  unfavorable. 

2.4  Experimental  Impact  Measurements 

The  calibrated  PMMA  bar,  was  used  in  subsequent  impact  measurements  using 
PMMA  and  PC  target  specimens.  Measurements  were  first  made  in  which  a  spherical 
pendulum  bob  served  as  the  projectile.  This  arrangement  paralleled  our  prior  work 
using  the  A1  bar,  and  ensured  good  control  of  the  impact  position,  and  the  accurate 
determination  of  the  impact  and  rebound  velocities.  These  velocities  are  needed  to 
compute  the  momentum  and  kinetic  energy  parameters.  Momentum  can  also  be 
determined  with  less  precision  from  the  recoil  swing  of  the  suspended  impact  bar. 

In  shifting  to  the  compressed  gas  gun,  the  projectile  velocity  was  determined  from 
the  time  interval  to  interrupt  two  light  beams  a  known  distance  apart.  The  projectile 
rebound  velocity  could  no  longer  be  measured  directly,  but  had  to  be  inferred  from 
the  recoil  swing  of  the  impact  bar.  For  this  work  a  second  strain  gage  was  added  the 
same  distance  from  the  first  as  the  first  is  from  the  specimen  end  of  the  bar.  The  two 
strain  gage  signals  were  recorded  at  a  common  time  base  at  0.5  /xs  intervals.  This 
second  gage  was  originally  intended  to  provide  a  further  check  on  the  validity  of  the 
previously  determined  transfer  function.  Ballistic  measurements,  using  the  gas  gun, 
were  made  in  duplicate  on  PC  and  PMMA  at  two  different  velocity  levels  using  the 

4.5  mm  ball  bearings  as  the  projectiles.  Very  good  reproducibility  was  experienced  be¬ 
tween  duplicate  runs. 

2.5  Analysis  of  Experimental  Results 

In  the  cases  of  impacts  using  glass  and  SiC  targets,  in  the  ideal  dispersionless  situa¬ 
tion  it  is  valid  to  invert  Eq.  1.1  to  yield 

x  -  u(t  -  t0)  -  ( P/K 0)^3  (Eq.  2.4) 

in  which 

P  -  AEt 
or 

t-  t0+l(AE/K0)W/UW 
-  t„  +  CtV3 

In  this  relationship  P  is  force,  A  is  the  bar  cross  sectional  area,  E  is  Young’s 
modulus,*  the  constant  in  Eq.  1.1,  in  the  impact  velocity,  C  the  proportionality 
constant  defined  by  the  equations,  e  the  observed  strain,  t  is  the  clock  reading,  and  t„ 

*  Henceforward  the  subscript  will  be  omitted  from  the  time-independent  Young's  modulus. 


is  the  cloctc  reading  at  the  time  of  initial  contact.  Thus,  plotting  c2/3  against  time  leads 
to  a  linear  relationship,  and  t'.ie  intercept  with  the  zero  strain  axis  fixes  t0.  Because  of 
the  lack  of  an  a  priori  theoretical  model  for  the  behavior  of  the  leading  edge,  and  be¬ 
cause  of  the  restrictions  on  the  functionality  of  the  pulse  shape  imposed  by  the 
transfer  function  approach,  our  working  hypothesis  was  that  the  toe  of  the  pulse 
should  increase  at  no  more  than  the  3/2  power  of  the  pulse  time.  Hence,  such  a  plot 
should  remain  useful  because  the  curve  at  t0  should  be  either  vertical  or  have  a  finite 
slope. 

However,  such  plots  of  the  new  data  obtained  using  the  gas  gun  did  not  indicate  a 
tendency  for  the  slope  to  approach  the  zero  strain  axis  with  a  finite  slope,  but  rather 
persisted  in  tending  towards  this  axis  with  zero  slope.  At  this  point  the  preconceived 
condition  of  Hertz-like  behavior  of  the  propagated  pulse  was  relaxed  and  instead  the 
data  was  futed  to: 

!-  t„+  CV/"  (Eq.  2.5) 

where  n  was  that  exponent  that  provides  the  best  empirical  fit  to  the  data,  and  C  is 
new  proportionality  constant.  A  good  fit  over  a  large  portion  of  the  start  of  the  pulse 
was  obtained  for  r.  =  4  for  the  data  from  the  first  gage  position  and  similarly  for  «  — 
6.3  from  the  second  gage  position.  This  result  suggested  that  more  of  the  start  of  the 
pulse  was  obscured  by  the  zero  signal  noise  than  had  been  realized.  These  plots  al¬ 
lowed  a  well-defined  extrapolation  to  the  zero  axis  within  experimental  scatter. 
Justification  that  these  power  functions  are  physically  significant  is  given  in  Sec¬ 
tion  5.4.1.  These  pulse  arrival  times  could  then  be  used  to  estimate  the  velocity  at 
which  the  leading  edge  advances.  This  result  was  found  to  be  very  close  to  the 
compressional  wave  velocity.  However,  most  of  the  pulse,  such  as  the  peak  portion 
travels  at  a  velocity  that  is  characteristic  of  the  extensional  wave  propagation  expected 
for  longitudinal  sound  propagation  in  a  cylindrical  rod. 

These  velocity  results  conflict  with  the  expectations  of  the  classical  theory  of  con¬ 
tinuous  wave  propagation  in  a  cylinder  discussed  in  1.1.2.  The  classical  theory  pre¬ 
cludes  sound  propagation  in  a  bar  in  excess  of  the  extensional  velocity.  This 
discrepancy  between  physical  expectation,  which  is  concordant  with  our  results,  and 
the  predictions  of  the  classical  continuous  wave  analysis  has  been  noted  by  Kolsky.(9) 

Furthermore,  if  the  correct  value  of  n  is  4  or  6.3,  then  there  is  no  mathematical 
possibility  of  finding  a  transfer  function  that  can  work  in  Eq.  2.2.  This  conclusion  car¬ 
ries  the  implication  that  the  Fourier  integral  app'oach,  using  the  above  classical  theory 
for  the  dispersion  correction  terms,  would  have  led  to  erroneous  results.  These  con¬ 
clusions  indicated  that  an  entirely  new  approach  to  the  dispersion  correction  problem 
was  required. 

2.6  Semi-Empirical  Approach  to  Pulse  Shape  Correction 

A  careful  examination  of  the  experimental  data  revealed  that  n  in  Eq.  2.5  and  prop¬ 
agation  velocities  of  the  leading  edge  and  the  pulse  maximum  were  constants.  The 
results  were  independent  of  the  material  impacted  or  the  impact  velocity.  The  ob¬ 
served  pulses  could  be  empirically  fitted  from  the  leading  edge  to  the  pulse  maximum 
with  little  error  by  an  equation  of  the  form: 

Y  —  Axn  —  Bxm  (Eq.  2.6) 

in  which  Fis  the  amplitude,  xis  the  time  since  the  first  arrival  of  the  pulse,  and  A ,  £, 
ft,  and  m  are  constants  to  be  determined.  The  requirement  that  the  peak  correspond 
in  magnitude  and  time  with  the  observed  value,  and  that  the  exponent  ft  at  the  leading 
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edge  also  match  that  found,  leaves  only  one  adjustable  constant  to  At  the  rest  of  the 
data.  The  position  of  the  peak  relative  to  the  leading  edge  is  predictable  front  the  two 
propagation  velocities.  The  dependence  of  n  on  position  can  be  inferred  in  pat  t  from 
the  measurements,  and  in  part  on  theoretical  grounds.  The  latter  is  based  on  the  fol¬ 
lowing  argument.  The  leading  edge  propagates  faster  than  the  rest  of  the  pulse.  This 
extension  can  be  accommodated  by  an  increase  in  the  power  dependence  on  t  with  in¬ 
creasing  distance  travelled.  The  shape  to  which  a  given  pulse  will  evolve  as  it  propa¬ 
gates  can  be  validly  determined  by  convolution.  If  the  transfer  function  appropriate  to 
this  process  leads  to  an  increase  in  the  power  dependence  with  increase  of  propagation 
distance,  then  conversely  the  starting  pulse  must  have  had  a  lesser  power  dependence. 
The  same  transfer  function  can  be  assumed  to  operate  over  any  portion  of  the  bar  of 
the  same  length.  Thus,  a  Axed  change  in  the  power  dependence  is  expected  over  a 
Axed  distance.  In  other  words  the  exponent  of  about  6.3  at  the  second  sensor  posi¬ 
tion,  and  of  about  4  at  the  Arst  sensor  position,  is  a  change  of  exponent  of  2.3,  This 
in  turn  implies  an  exponent  of  about  1.7  at  the  bar  end  and  1.5  at  the  impact  site. 
(This  is  the  value  computed  using  the  Hertz  elastic  model.)  Thus,  the  essence  of  the 
transfer  function  concept  has  been  retained,  except  that  convolution  as  a  mathematical 
procedure  has  been  supplemented  by  a  more  qualitative  procedure. 

This  overall  approach  to  pulse  correction  is  discussed  in  depth  in  Section  5.4.  By 
this  means  the  observed  pulses  were  corrected  to  infer  the  shape  of  the  initial  impact 
pulse.  The  corrected  pulse  was  then  introduced  into  the  procedures  for  determining 
the  penetration  behavior. 

2.7  Derived  Material  Response  Behavior 

The  corrected  pulse  data  were  analyzed  to  deduce  the  force-penetration  depth, 
penetration  depth-time,  and  local  energy-penetration  deoth  dependences.  The  results 
on  the  PMMA  bar  using  the  pendulum  impacts  largely  conArmed  the  results  previous¬ 
ly  obtained  using  the  A1  bar.  The  new  procedure  appears  to  have  eliminated  the  ar¬ 
tifacts  noted  in  our  previous  studies  in  the  high  velocity  impacts.  The  results  on 
PMMA  using  the  compressed  gas  gun  appear  to  be  extensions  of  the  data  produced 
using  the  mechanical  test  machine,  whereas  those  results  on  PC  indicate  an  anomalous 
softening  with  increasing  impact  velocity.  The  results  are  presented  in  Section  6.1  and 
6.2.  Both  materials  exhibit  evidence  of  load  relaxation  during  the  course  of  the  pene¬ 
tration.  The  peak  loads  occur  prior  to  reaching  the  peak  penetration  depths. 

2.8  Modelling  of  the  Material  Behavior 

A  visco-elastic  model  of  the  material  behavior  was  developed  which  incorporates 
the  elements  of  the  model  used  to  explain  velocity-dependent  relaxation  behavior 
given  by  Eq.  1.5.  The  model  is  aimed  at  accounting  for  the  observed  penetration- 
time-force  inter-relationships,  at  a  better  understanding  of  how  the  molecular  and 
structural  responses  are  influenced  by  the  initial  impact  conditions,  and  at  trying  to  ac¬ 
count  for  the  differences  between  the  polymers.  The  model  is  developed  in  terms  of  a 
system  of  elastic  responses  and  velocity-time-dependent  viscosity  responses.  These 
responses  act  partly  in  parallel  and  partly  in  series.  In  PC  it  appears  that  the  initial 
response  is  almost  a  purely  parallel  behavior.  In  PMMA  the  behavior  is  about  25% 
parallel  and  75%  series. 


r-’M-T, rtldfafr  ’"MM 
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3.  MATERIALS  AND  SPECIMEN  PREPARATION 

The  PC  and  PMMA  materials  were  unchanged  from  those  used  in  the  previous 
study.  Specimens  were  used  as  cut  out  of  sheet  stock,  making  sure  that  all  burrs  were 
removed.  The  standard  specimen  size  was  25  mm  O.D.  x  12.5  mm  The  other  target 
materials  are  described  below. 

3.1  Polymethylmethacrylate 

Specimens  were  cut  from  a  single  sheet  of  Type  G  Plexiglass R  PMMA  (produced 
by  Rohm  and  Haas).  Continuing  with  our  former  practice,  the  sheet  was  marked  off 
in  squares  for  cu  ing  and  a  record  kept  so  that  the  original  location  of  each  specimen 
in  the  sheet  could  be  identified. 

3.2  Polycarbonate 

The  PC  resin  was  manufactured  by  the  General  Electric  Company  and  is  designated 
as  LexanR  resin  general  purpose  glazing  sheet,  Type  9034-1 12.  A  single  sheet  of 
12.5  mm  thick  material  was  used,  and  specimens  marked,  and  records  kept  as  for  the 
case  of  PMMA. 

3.3  Inorganic  Glass 

A  sheet  of  Pyrex*  borosilicate  glass  (made  by  Corning  Glass  Works)  12.5  mm 
thick  was  ground  into  125  mm  diameter  disks  for  use  in  calibrating  the  SiC>2  impact 
bar. 

3.4  Silicon  Carbide 

A  piece  of  fully  dense,  hot  pressed  silicon  carbide  produced  by  the  Norton  Co.  was 
ground  into  a  cylinder  having  a  12.5  mm  diameter  and  a  11.9  mm  height.  This  height 
corresponds  to  a  calculated  acoustic  reverberation  time  of  2  fis. 
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4.  IMPACT  BAR  STUDIES 

The  experimental  arrangement  of  the  impact  bar  is  in  many  respects  unchanged 
from  the  previous  contract  study.  The  main  differences  are  that  two  new  bars  were 
made,  one  of  vitreous  silica,  the  other  of  PMMA  having  diameters  of  12.5  mm  and 
25  mm,  respectively.  The  wire  strain  gages  were  replaced  by  semiconductor  gages. 
The  resultant  improvement  in  signal  strength  eliminated  the  need  for  pre-amplification 
of  the  signal  ahead  of  the  Nicolet  Model  2090-/  storage  oscilloscope. 

The  specimens  were  mounted  on  one  end  of  the  180  cm  long  instrumented  impact 
bar.  The  specimen  diameters  either  matched  or  were  less  than  the  bar  diameter  to 
avoid  reflections  at  the  overhang.  The  arrangement  is  shown  schematically  in  Fig.  4A. 

Strain  gages  on  the  bar  connected  to  a  Wheatstone  bridge  and  the  output  to  a  digi¬ 
tal  storage  oscilloscope  enabled  the  strain  pulse  produced  by  the  impact  to  he  recorded. 
Mathematical  processing  of  the  pulse  data,  using  a  simplified  one-dimensional  model 
for  the  stress  distribution  in  the  impact  bar,  allows  the  force,  penetration  depth,  striker 
velocity,  various  energy  terms  and  the  momentum  transfer  to  be  computed.  This  pro¬ 
cedure  is  given  in  the  Appendix.  The  polymer  PMMA  was  the  final  selection  as  the 
impact  bar  material  because  of  its  very  good  acoustic  impedance  match  with  the  target 
specimens.  Secondly,  the  small  value  of  its  Young's  modulus  resulted  in  a  large  strain 
for  a  given  force.  This  was  favorable  with  respect  to  signal-to-noise  considerations. 

The  bar  was  supported  at  two  places  along  its  length  by  two  bifilar  fine  wires  in  the 
general  shape  of  a  V.  This  allows  the  bar  to  swing  freely  along  its  axis,  while  substan¬ 
tially  preventing  sidewise  displacement.  A  similar  arrangement  was  used  to  support 
the  striker.  From  the  amplitude  of  the  swing  of  the  bar  after  the  specimen  is  impacted 
by  the  striker  and  from  the  initial  position  of  the  striker  before  its  release,  the  mo¬ 
menta  and  kinetic  energies  can  be  computed.  Two  striker  geometries  were  used.  One 
was  a  19  mm,  28.3  g  steel  hall.  The  other  was  a  42  g  steel  cylinder  with  one  end  in 
the  form  of  a  nose  terminating  in  a  4.5  mm  diameter  spherical  cap.  The  latter 
geometry  provided  continuity  with  prior  indentation  work  and  with  the  air  gun  work 
using  4.5  mm  projectiles.  This  pendulum  striker  arrangement  was  convenient  for  im¬ 
pact  velocities  up  to  about  2.5  m/s. 

Higher  impact  velocities  were  produced  by  firing  a  4.5  mm  steel  ball  from  a 
modified  air  rifle  using  a  controlled  pressure  of  He  as  the  propellant.  The  gun  was 
mounted  on  an  I-beam  equipped  with  aiming  adjustments  as  well  a?  with  a  pair  of 
photocells  placed  a  known  distance  apart.  The  velocity  of  the  projectile  was  deter¬ 
mined  on  an  interval  timer  triggered  by  the  interruption  of  the  light  beams  aimed  at 
the  photocells.  Because  of  degradation  of  the  accuracy  of  the  trajectories  at  both  high 
and  low  propellant  pressures,  the  impact  velocities  were  restricted  to  between  120  and 
220  m/s.  The  projectile  was  aligned  to  be  coincident  with  the  axis  of  the  bar  and  the 
center  of  the  target  specimen  mounted  on  the  suspended  impact  bar. 

In  the  case  of  the  gas  gun  experiments  two  sets  of  strain  gages  were  mounted  on 
the  PMMA  impact  bar,  one  set  15  cm  from  the  end  of  the  bar,  and  the  second  set 
30  cm  from  the  end.  The  pulse  signals  detected  by  each  of  the  two  sets  could  be 
recorded  simultaneously  in  separate  channels  of  the  oscilloscope. 

4.1  Impact  Measurements  Using  Vitreous  Silica  Bar 

The  vitreous  silica  bar  was  used  at  the  beginning  of  this  study  because  its  unusually 
small  Poisson  ratio  results  in  an  advantageously  small  dispersion.  Impact  measure¬ 
ments  were  made  on  Pyrex  glass  and  hot-pressed  silicon  carbide  specimens  which  were 
impacted  by  a  19  mm  steel  ball  bearing  at  a  velocity  of  1.77  m/s.  Preliminary  results 
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made  on  specimens  having  a  larger  cross  section  than  the  bar,  when  compared  with 
later  results  showed  the  importance  of  using  targets  that  match  the  cross  section  of  the 
bar.  This  is  shown  in  Figs.  4.1  A  and  B. 

4.2  Use  of  Vitreous  Silica  as  a  Striker 

The  small  degree  of  dispersion  in  the  vitreous  silica  bar  made  it  advantageous  for 
use  as  a  means  of  calibrating  the  PMMA  bar.  For  the  reasons  discussed  in  Sec¬ 
tion  2.3  SiC2  was  unsuitable  for  use  as  the  bar  material  in  combination  with  polymer 
targets  in  high  velocity  impact  experiments.  However,  by  mounting  a  glass  hemi¬ 
sphere  on  the  end  of  the  Si02  bar,  it  could  be  made  into  a  particularly  ideal  striker. 
The  strain  gage  on  the  Si02  bar  would  sense  a  relatively  distortion-free  signal  arising 
from  its  impact  with  the  target  on  the  PMMA  bar.  Hence,  the  pulse  as  detected  in  the 
Si02  bar  could  be  directly  compared  with  that  detected  in  the  more  dispersive  PMMA 
bar.  The  nominal  similarity  between  the  two  signals  of  the  same  impact  event  are 
seen  in  Fig.  4.2A.  Note  that  the  times  do  not  coincide  because  of  the  differences  in 
the  propagation  velocities  and  in  the  positions  of  the  gages  in  the  two  bars.  Note  also 
the  more  gradual  onset  of  the  leading  edge  signal  in  the  case  of  the  PMMA  bar. 

4.3  Impact  Measurements  Using  the  PMMA  Bar 

Using  the  42.0S  g  striker  that  is  fitted  with  a  4.5  mm  diameter  tip,  impact  measure¬ 
ments  were  made  in  duplicate  on  PC  and  PMMA  specimen  targets.  The  impact  veloc¬ 
ity  was  about  2.1  m/s  in  all  cases.  A  representative  example  is  given  in  Fig.  4.3A. 
There  was  good  agreement  between  duplicate  measurements  as  shown  in  the  table. 
The  pulses  were  similar  in  behavior  and  appeared  to  be  nearly  classical  in  form.  The 
pulse  rose  to  a  maximum  and  diminished  again  without  crossing  the  zero  strain  axis. 
The  only  slight  irregularity  is  that  the  strain  level  remained  positive  after  the  pulse  had 
apparently  passed.  This  may  be  due  to  a  slight  creep  of  the  PMMA  impact  bar  under 
the  pulse  load. 

A  second  group  of  experiments  was  made  using  the  gas  gun  to  produce  the  impacts 
on  the  PC  and  PMMA  targets.  Duplicate  experiments  were  run  at  nominally  125  and 
220  m/s  on  both  PC  and  PMMA  targets.  The  pulses  at  the  two  gage  positions  were 
recorded  in  each  case.  The  duplicate  runs  gave  remarkably  similar  signals,  even  in¬ 
cluding  the  complex  structure  of  the  trailing  portion  of  overall  pulse.  This  can  be  seen 
by  comparing  Fig.  4.3  B  with  C.  The  pulse  received  at  the  second  gage  typically  had 
less  fine  structure  and  displayed  larger,  better  developed  lower  frequency  oscillations 
as  may  be  noted  by  comparing  Fig.  4.3  C  with  D.  While  the  pulse  pattern  for  PC  and 
PMMA  were  different  in  detail,  they  were  comparable  in  general  character.  This  can 
be  verified  by  comparing  Fig.  4.3  B  with  E.  The  most  noticeable  difference  is  the 
somewhat  longer  period  for  the  initial  half  wave  in  the  case  of  PC  relative  to  PMMA. 
Finally,  there  appeared  to  be  little  general  effect  resulting  from  increasing  the  impact 
velocity  from  125  m/s  to  220  m/s  other  than  increasing  the  amplitude  of  the  pulse. 
This  is  shown  in  Fig.  4.3  F  relative  to  B. 


S.  CORRECTION  OF  IMPACT  PULSE  DATA  FOR  DISPERSION 

This  is  the  most  crucial  portion  of  the  work  because  of  its  effect  on  the  calculations 
that  reveal  the  impact  response  of  the  target  specimen  material.  If  the  force-time  rela¬ 
tionship  at  the  site  of  the  impact  were  known  exactly  ,  then  the  force  could  be  used  to 
compute  the  projectile  deceleration  and  the  force-depth  relationship  could  be  deter¬ 
mined  exactly.  Thus,  the  problem  is  to  transform  the  pulse  data  observed  at  a  later 
time  and  at  a  position  remote  from  the  impact  site,  into  the  original  pulse  shape.  For 
these  purposes  the  correct  assignment  of  the  start  of  the  pulse,  and  the  correct  quanti¬ 
tative  description  of  the  force  during  the  initial  stages  of  the  projectile  penetration  are 
vital  to  determining  the  characteristic  visco-elastic  properties  of  the  target  material. 
The  correct  shape  of  the  entire  pulse  is  needed  to  reveal  the  relaxation,  the  maximum 
penetration  and  force,  the  energy  absorption  processes  as  well  as  the  rebound 
response. 

The  correction  procedure  is  first  to  concentrate  on  the  leading  edge  behavior  to 
define  its  arrival  time  and  shape.  Then  the  overall  pulse  shape  and  the  rate  at  which 
the  various  regions  of  the  pulse  propagate  are  considered.  These  differences  in  the 
propagation  velocity  are  responsible  for  the  pulse  distortion.  The  quantitative  descrip¬ 
tion  of  the  evolution  of  the  shape  is  then  inverted  to  deduce  the  shape  at  a  preceding 
time  and  location,  specifically  the  time  and  place  of  the  impact.  Second  order 
differences  determined  earlier  in  the  process  are  reintroduced  as  fine  adjustments. 
This  completes  the  correction  of  the  observed  pulses. 

5.1  Detection  and  Characterization  of  the  Leading  Edge 

Leading  edge  arrival  could  be  observed  at  three  independent  sensor  locations  viz. 
the  Si02  bar  gage,  and  the  two  gages  on  the  PMMA  bar.  However,  because  the 
dispersion  in  the  Si02  bar  was  so  small,  the  start  of  the  pulse  was  virtually  directly  ob¬ 
servable  in  the  special  case  of  the  collision  between  these  two  bars. 

5. 1. 1  Vitreous  Silica  Bars 

In  the  impact  of  the  191mm  steel  ball  with  the  Pyrex  specimen  mounted  on  the 
Si02  impact  bar  the  conditions  for  an  elastic  impact  were  satisfied.  Hence,  as  previ¬ 
ously  shown  in  Eq.  2.4  the  (2/3)  power  of  the  strain  should  increase  linearly  with  time 
starting  at  the  instant  of  first  contact.  The  plot  in  Fig.  5.1A  for  the  experimental 
results  of  the  average  of  the  two  duplicate  runs  show  a  slight  departure  near  the  very 
start  of  the  impact.  The  deviation  between  the  extrapolated  impact  time  and  the  prob¬ 
able  actual  time  is  about  1  ns.  This  indicates  that  the  signal  from  this  bar  can  be  used 
without  correction  to  an  accuracy  of  about  1  ns  for  pendulum  impacts  such  as  those 
used  in  the  present  experiments. 

5.1.2  PMMA  Impact  Bar 

5. 1.2.1  Vitreous  Silica  Bar  Collision  with  PMMA  Bar 

With  the  knowledge  that  the  signal  from  the  Si02  impact  bar  is  a  good  representa¬ 
tion  of  the  actual  impact  process,  it  is  interesting  to  examine  the  leading  edge  behavior 
for  this  experiment.  A  representative  plot  of  the  2/3  power  of  the  observed  strain  vs 
time  in  the  PMMA  bar  is  shown  in  Fig.  5.1  B.  Whereas  the  signal  in  the  Si02  is  seen 
to  be  linear  within  experimental  error,  the  corresponding  curve  for  the  PMMA  signal 
which  is  linear  some  distance  above  the  zero  strain  axis,  clearly  does  not  approach  the 
abscissa  linearly.  The  start  of  the  pulse  is  difficult  to  define,  but  cursory  examination 
suggests  that  it  may  occur  at  least  10  ps  prior  to  the  time  suggested  by  extrapolating 
the  more  linear  portion  back  to  the  abscissa. 


However,  if  the  data  is  represented  instead  by  the  empirical  Eq.  2.5  then  the  linear 
behavior  shown  in  Fig.  5.1  C  results  for  n  *  4.  A  duplicate  run,  not  shown,  for  the 
two  bar  impact  gives  the  same  linear  behavior  for  this  value  of  n  There  is  increase  in 
scatter  near  the  zero  axis  due  to  the  magnification  that  occurs  in  small  quantities,  such 
as  random  noise,  when  their  4th  root  is  extracted.  Values  of  the  4th  root  less  than 
about  0.08  are  attributable  to  noise,  in  spite  of  the  long  extrapolation  to  the  zero 
strain  axis,  the  intercept  can  be  defined  to  ±  1  fis.  If  this  extrapolation  is  indeed  the 
true  arrival  time  of  the  pulse,  then  much  of  the  start  of  the  pulse  is  obscured  by  noise 
before  it  can  be  unambiguously  detected. 

It  should  be  noted  that  the  non-linear  portion  of  the  2/3  power  piot  as  well  as 
much  of  the  linear  portion  is  linearized  by  the  4th  root  plot.  The  band  of  noise  sur¬ 
rounding  the  zero  strain  axis,  made  apparent  by  the  latter  plot,  is  not  as  evident  in 
the  2/3  power  plot. 

5.1. 2.2  Pendulum  Impacts 

The  start  of  the  pulses  for  the  pendulum  impacts  of  PMMA  and  PC  specimens  by 
the  smaller  4.5  mm  diameter,  42  g  striker  has  been  similarly  analyzed.  The  pulses 
sensed  at  the  gage  15  cm  from  the  bar  end  (16.2  cm  from  the  location  of  the  impact), 
in  all  cases  show  a  non-linear  behavior  at  small  strains  when  examined  by  the  2/3 
power  plot.  As  before,  identification  of  the  start  of  the  pulse  is  difficult  as  shown  in 
Figs.  5.1  D-E.  However,  when  the  4th  power  of  strain  is  plotted  instead,  an  approxi¬ 
mately  linear  behavior  is  again  seen  in  Figs.  5.1  F-G.  The  increasing  prominence  of 
noise  at  small  strain  levels  is  again  evident. 

5. 1.2.3  Gas  Gun  Projectile  Impacts 

The  leading  edge  of  the  pulse  received  at  the  15  cm  gage  position  for  the  impacts  at 
120  and  220  m/s,  were  similar  to  those  from  the  pendulum  impacts.  Representative 
plots  of  the  4th  root  of  strain  are  given  in  Figs.  5.1  H-I.  However,  the  leading  edge 
signals  detected  at  the  second  gage,  located  30  cm  from  the  end  of  the  bar  could  only 
be  linearized  if  6th  root  of  the  strain  was  used  instead.  The  appropriate  root  was 
found  by  varying  the  exponent  n  one  unit  at  a  time  until  a  linear  curve  was  obtained. 
These  results  in  Figs.  5.1  J-K  also  show  the  relative  prominence  of  noise  is  increased 
even  more.  This  increases  the  span  of  the  extrapolation  needed  to  define  the  start  of 
the  pulse  signal.  Consequently,  there  is  greater  uncertainty  in  the  extrapolated  value. 
A  procedure  for  refining  the  estimate  of  the  exponent  is  given  in  5.3. 

5.2  Propagation  Velocity 

In  the  experiments  involving  the  two  bar  collision  and  in  the  high  velocity  impact 
experiments  the  pulse  shape  at  two  bar  locations  was  simultaneously  recorded  on  a 
common  time  base.  Thus,  the  propagation  velocity  can  be  determined  for  the  leading 
edge,  the  pulse  maximum,  the  time  of  the  zero  crossing,  and  for  other  definable  por¬ 
tions  of  the  impact  pulses.  In  the  case  of  the  Si02  bar  the  extensional  and  compres- 
sional  wave  velocities  are  almost  identical  being  5,760  and  5,968  m/s,  respectively. 
The  gage  on  the  Si02  bar  was  7.5  cm  from  the  bar  end  (8.7  cm  from  the  position  of 
the  specimen  impact).  Therefore  the  propagation  time  from  the  impact  site  to  the 
sensor  is  about  15.0  ±  0.5  #xs.  Since  there  was  a  common  time  base  for  the  signals  in 
the  two  bars,  subtracting  15  /xs  from  the  clock  reading  when  the  pulse  first  arrived  at 
the  gage  in  the  Si02  bar  represents  the  clock  reading  at  the  start  of  the  impact.  The 
4th  root  plot,  discussed  above,  was  then  used  to  estimate  the  arrival  time  of  the  pulse 
in  the  PMMA  bar  at  the  gage  16.2  cm  from  the  impact  site.  From  this  data  the  propa¬ 
gation  velocity  in  the  PMMA  bar  for  the  hading  edge  was  calculated  to  be  2660 
±  30  fi/s.  The  error  represents  the  variance  between  two  duplicate  runs.  The  very 


19 


long  duration  of  the  two  bar  impact  pulse  precluded  making  velocity  measurements 
on  other  parts  of  the  pulse. 

A  more  direct  determination  of  the  propagation  velocities  is  possible  when  the  two 
sensors  are  located  on  the  same  bar  and  there  is  again  a  common  time  base.  In  the 
high  velocity  impact  experiments  the  pulse  arrival  times  at  the  first  gage  position  were 
estimated  on  the  basis  of  the  4th  root  plots,  while  those  at  the  second  position  15  cm 
from  the  first  were  estimated  on  the  basis  of  the  6th  root  plots.  The  procedure  for 
determining  the  velocities  is  shown  schematically  in  Fig.  5.2  A.  The  results  are  given 
in  Table  5A.  Note  that  the  velocity  results  depend  mainly  on  the  position  within  the 
pulse  rather  than  on  how  the  pulse  was  initially  produced. 


Table  SA 

CHARACTERISTICS  OF  PROPAGATION  OF  IMPACT  PULSE  BETWEEN 
15  cm  AND  30  cm  SENSOR  LOCATIONS  ON  2.5  cm  DIAMETER  PMMA  BAR 


Pulse 

Location 

Impact 

Velocity 

m/s 

Propagation  Time  Between  Gages 

US 

Mean 

Velocity 

PC  Target 

PMMA  Target 

m/s 

Leading  edge 

120 

56.85±.25 

54.75±.7S 

2710±  80 

220 

54.50±2.4 

55.55±.95 

Half  Max.  Ht. 

120 

72.45  ±.05 

72.60±.00 

2070±  10 

(Ascending) 

230 

72.40±1 

72.70±.25 

Maximum 

110 

75.05±.05 

74.90±.05 

2000±5 

220 

75.10±.l 

75. 00 ±.05 

Half  Max.  Ht. 

110 

75.95  .05 

77.12±.03 

1960±15 

(Descending) 

220 

75.88 

76.90±.35 

Zero  Crossing 

110 

78.50±1.2 

77.75±.05 

1920±  25 

220 

78.00±0.7 

77.80±.05 

5.3  Analytical  Representation  of  Overall  Pulse 

The  total  pulse,  including  the  complex  oscillations  of  the  tail  portion  has  a  very 
complicated  shape.  These  oscillations  are  due  to  internal  reflections,  dispersion,  the 
excitation  of  propagation  modes  other  than  the  simple  longitudinal  one,  and  to  the 
later  arrival  of  shear  and  surface  waves.  Axial  momentum  p  is  calculable  from  the  in¬ 
tegral, 

pit)  -  f  Pdt  (Eq.  5.1) 

-  EA  f  fdt 
*'o 

where  Fis  force,  £is  the  Young’s  modulus  of  the  bar,  A  its  cross  section,  e  the  axial 
strain  in  the  bar,  and  t  is  time.  When  the  integration  of  the  measured  strain  over  time 
was  performed,  it  was  found  that  the  oscillations  beyond  the  first  zero  crossing  very 
nearly  cancel  and,  therefore,  are  not  involved  in  the  net  axial  momentum  transfer. 
Hence,  we  are  effectively  concerned  only  with  the  pulse  shape  up  to  the  time  that  the 
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strain  first  crosses  the  abscissa.  This  resuits  in  an  enormous  simplification  of  the  anal¬ 
ysis.  In  any  case,  it  is  this  first  portion  of  the  pulse  that  contains  most  of  the  informa¬ 
tion  about  the  material  behavior  and  the  impact  event. 

As  discussed  in  5.1  and  5.2  the  arrival  time  of  the  leading  edge  can  be  “defined" 
once  the  appropriate  n  is  found  that  results  in  a  linear  curve.  (The  argument  for  the 
validity  of  this  procedure  for  defining  the  true  shape  and  pulse  arrival  time  is  given  in 
5.4.)  Thus,  the  start  of  the  pulse  is  given  by: 

strain  -  constant  •  tn 

where  r  is  the  time  since  first  arrival  of  the  pulse.  Other  characteristics  of  the  pulse 
are  its  maximum  strain  e(max)  and  the  elapsed  time  r(max)  at  this  point.  For  con¬ 
venience  we  define 

S  -  t/  t(max)  (Eq.  5.2a) 

and 

y  -  e/e  (max)  (Eq.  5.2b) 

Then  the  following  equation  has  the  characteristics  of  having  the  proper  dependence 
of  strain  on  time  at  the  start  of  the  pulse,  and  the  proper  value  and  location  of  the 
maximum 


y  - 


sn-  sm 

n  —  m 


(Eq.  5.2c) 


in  which  m  is  an  adjustable  constant,  determined  by  ar.  error  minimization  in  compari¬ 
son  with  the  actual  data.  In  fact  using  this  equation,  which  is  the  normalized  form  of 
Eq.  2.6,  both  n  and  m  can  be  varied  to  give  the  best  fit.  This  results  in  a  corrected  es¬ 
timate  for  n. 

The  preceding  equation  models  the  high  velocity  pulse  data  up  to  the  maximum  to 
within  an  error  of  a  few  per  cent  absolute,  as  shown  in  Fig.  5.3  A,  corresponding  to 
about  10%  relative.  Because  of  the  small  variation  in  propagation  velocity  noted  for 
the  descending  portion  of  the  pulse,  there  is  only  a  small  change  in  shape  as  the  pulse 
travels.  Therefore,  this  portion  of  the  pulse  has  not  been  modeled.  Rather  the  time 
scale  is  adjusted  according  to  the  pulse  position  and  bar  location. 

5.4  Pulse  Correction  Methodology 

This  section  is  closely  related  to  2.1  and  2.6,  which  are  recommended  as  back¬ 
ground.  Many  of  the  ideas  and  basic  information  were  shown  schematically  in 
Fig.  5.2  A. 

5.4.1  Conceptual  Basis 

Introduction  of  a  variable  exponent  into  the  empirical  description  of  the  toe  region 
of  an  impact  pulse,  led  to  the  interesting  result  that  at  a  given  sensor  position  the 
same  exponent  was  applicable  to  a  wide  range  of  pulses.  The  exponent  was  about  4  at 
the  first  gage  position,  and  somewhat  greater  than  6  at  the  second  gage  position. 
These  exponents  led  to  consistent,  well  defined  values  for  the  arrival  times  of  the 
pulses.  The  value  for  the  propagation  velocity  of  the  leading  edge  based  on  the  data 
from  the  two  gage  positions  on  the  PMMA  bar  was  about  2610  m/s.  This  result  com¬ 
pares  closely  with  the  value  2660  m/s  found  for  the  two  bar  impacts  using  the  data 
from  the  gage  on  the  Si02  bar  and  that  on  the  PMMA  bar.  These  two  velocities 
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values  are  remarkably  close  to  2680  m/s  which  is  the  reported  longitudinal  sonic  ve- 
locity  for  PMMA.  (The  longitudinal  velocity  represents  greatest  sonic  propagation  ve¬ 
locity  possible.)  The  velocities  found  for  the  pulse  peak  and  the  descending  portion 
more  closely  approach  the  reported  extensions!  velocity  1840  m/s.  The  internal  con¬ 
sistency  of  this  body  of  data  is  supportive  of  the  hypothesis  of  an  underlying  physical 
basis.  The  reason  that  the  exponents  n  should  depend  on  position  in  the  way  noted 
follows  from  the  following  plausibility  argument. 

We  assume  that  the  effects  of  dispersion,  whatever  their  origin,  are  cumulative  and 
that  all  parts  of  the  impact  bar  are  equivalent  in  their  dispersion  characteristics,  and 
that  the  effect  of  dispersion  can  be  determined  by  convolution.  Recalling  Eq.  2.1,  if 
Fit)  is  the  original  signal,  and  Git)  is  its  form  after  travelling  in  a  dispersive  medium, 
then 

Git)  -  f  Fir)Hit-r)dr  (Eq.  2.1) 

**o 

where  H  is  a  function  that  contains  the  dispersion  characteristics. 

Consider  the  asymptotic  representation  of  ^the  pulse  at  small  values  of  elapsed  (local) 
time  to  be  given  by  / 

F\$(t)  -  tf,5r4  (Eq.  5.3) 

and 

F30it)~  K30t*3  (Eq.  5.4) 

where  the  subscripts  refer  to  the  sensor  location  in  cm  and  A!  is  a  constant.  Thus,  in 
the  sense  just  discussed  for  purposes  of  Eq.  2.1 

Fit)  -  Fls(t)  (Eq.  5.5) 

and 

Git)  -  F30(t)  (Eq.  5.6) 

Dropping  the  proportionality  constants  one  can  write 

F30it)  -  Fxsit)t2>  (Eq.  5.7) 

This  result  would  have  resulted  from  convolution  provided  Hit )  proportional  to  /u, 
at  least  for  small  t  By  inspection  one  can  see  that  the  effect  of  Hit)  is  to  multiply 
Fit)  by  r2-3  when  the  sensor  location  is  moved  15  cm  away  from  the  source.  This  in¬ 
dicates  one  can  write  generally 

Git)  -  Fit)  t 23  (Eq.  5.8) 

for  any  such  shift.  Conversely  this  can  be  written 

Fit)  -  Git)  t  2  3  (Eq.  5.9) 

which  states  that  the  signal  15  cm  towards  the  source  has  the  functional  behavior  of 
the  local  pulse  divided  by  t2-3.  Thus,  if  the  signal  at  the  first  gage  is  given  by 
Eq.  (5.3),  then  the  signal  at  the  end  of  the  bar  is 

F0it)  -  Kxst4f23  -  (Eq.  5.10) 
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However,  the  impacted  surface  of  the  specimen  is  1.2  cm  from  the  end  of  the  bar  so 
that  the  exponent  at  the  site  of  the  impact  would  be  expected  to  be  less  than  1.7.  It 
was  probably  not  a  coincidence  that  in  the  two  bar  collision  experiment,  in  which  the 
power  dependence  on  fof  the  initial  pulse  was  1.5,  the  pulse  observed  at  the  first  gage 
had  an  exponent  of  4.  The  Hertz  elastic  model  Eq.  1.1  which  predicts  a  1.5  exponent 
dependence  on  penetration  depth,  is  expected  to  remain  approximately  valid  for 
visco-elastic  materials.  Therefore,  a  constant  value  of  4  for  the  exponent  n  at  the 
15  cm  gage  location,  and  of  6.3  at  the  30  cm  location  is  not  only  predicted  but 
demanded  by  convolution.  This  completes  the  arguments  in  support  of  the  universali¬ 
ty,  constancy  and  validity  of  the  empirically  noted  pulse  behavior. 

Although  the  convolution  procedure  is  valid  for  computing  the  subsequent  pulse 
shape  (because  in  this  case  the  process  leads  in  a  higher  power  dependence  on  /),  the 
inverse  is  mathematically  undefined.  Nevertheless,  Eq.  5.9  does  provide  a  rationale 
for  deducing  at  least  what  the  leading  edge  behavior  must  have  been  at  some  prior 
time  and  location.  We  shall  use  the  principle  contained  in  Eq.  5.9  as  one  of  the  foun¬ 
dation  elements  in  transporting  a  pulse  backwards  in  time,  i.e.  to  correct  it  for  disper¬ 
sion. 

5. 4.2  Operational  Procedure 

The  primary  assumptions  for  reconstructing  the  original  pulse  from  the  observed 
impact  pulse  data  are  (1)  that  the  power  law  description  of  Eq.  5.2c  remains  valid  for 
describing  the  shape  of  the  pulse  up  to  its  maximum,  (2)  that  the  propagation  veloci¬ 
ties  determine  the  time  shift  of  the  position  of  the  peak  relative  to  the  start  of  the 
pulse,  (3)  that  the  change  in  the  shape  of  the  leading  edge  is  given  by  the  prescribed 
change  in  the  power  dependence  on  t,  and  (4)  that  conservation  of  momentum 
defines  the  amplitude  scale  of  the  initial  pulse. 

Because  the  dependence  of  the  leading  edge  shape  on  t  increases  by  a  power  of  2.3 
for  a  15  cm  length  of  bar  travel  it  follows  that  the  power  increase  is  proportional  to 
distance  at  the  rate  j8  =  2.3/15  per  cm.  Thus,  if  the  observed  initial  power  depen¬ 
dence  at  the  first  sensor  is  n}5  then  the  assumed  dependence  at  the  site  of  the  impact, 

16.2  cm  forward  of  the  sensor  is 

»o  -  "is  “  0  •  16.2  (Eq.  5.11) 

The  difference  between  the  propagation  velocity  at  the  leading  edge  v0  and  the 
pulse  maximum  vm  leads  to  a  progressive  stretching  of  the  pulse.  Defining  the  period 
A(x)  of  the  pulse  as  the  interval  from  the  start  to  the  peak  at  a  location  xcm  from  the 
impact  site,  then  at  the  first  sensor  location  the  period  is  given  by: 

A  (16.2)  -  A(0)  +  16.2  •  (-* - -)  (Eq.  5.12a) 

V*  V0 

or 

A(0)  -  A (16.2)  -  16.2  (J - l-)  (Eq.  5.12b) 

Vm  V0 

This  new  value  of  the  exponent  n0  from  Eq.  5.11  is  inserted  into  Eq.  5.2c,  in  place 
of  n  The  corresponding  value  of  m  is  determined  by  the  two  requirements  that  at  the 
pulse  maximum,  i.e.  S  =  1,  the  slope  must  not  only  be  zero  but  the  curvature 
(second  derivative)  must  be  continuous  with  that  on  the  descending  part  of  the  pulse. 
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This  leads  to  the  relationship 


w  +  m  -  1  _  n0+  m0- 
A(16.2)2  A  (0)2 


(Eq.  5.13) 


which  in  turn  defines  the  value  of  m0  to  be  introduced  into  Eq.  5.2c.  This  result  as¬ 
sumes  that  the  velocity  over  the  descending  slope  is  nearly  constant  over  its  span,  so 
that  the  change  in  shape  was  neglected. 


In  synthesizing  the  calculated  impact  pulse  at  the  specimen  surface,  the  time  scale 
is  reintroduced  into  Eq.  5.2c  by  use  of 

r-S  A(O)  (Eq.  5.14) 

The  descending  portion  is  then  joined  to  the  ascending  portion  by  assuming  that  the 
velocity  varies  linearly  between  that  at  the  peak  and  that  the  zero  crossing  Vcations. 
The  error,  such  as  shown  in  Fig.  5.3A  between  the  modelled  portion  of  ihe  ;<se,  and 
the  experimental  data  is  added  back  onto  the  translated  pulse  at  the  new  time  positions 
determined  by  Eq.  5.14.  In  this  way  the  “fine  structure”  of  the  experimental  data  is 
preserved.  This  defines  the  corrected  pulse  shape  in  terms  of  amplitudes  relative  to 
the  peak  amplitude.  The  absolute  amplitude  is  fixed  by  introducing  a  multiplier  into 
Eq.  5.1  to  satisfy  the  condition  of  conservation  of  momentum. 

The  procedure  as  outlined  is  equally  applicable  to  any  shift  in  the  pulse  detection 
site.  Thus,  the  pulse  data  obtained  at  the  second  sensor  location  can  be  used  to  com¬ 
pute  the  pulse  shape  expected  at  the  first  sensor. 

Two  computer  programs  were  written  to  perform  the  pulse  correction.  The  first 
one  was  to  provide  an  accurate  empirical  description  of  the  observed  pulse.  The 
second  uses  that  description  to  “translate”  the  pulse  forward  in  time  and  position  to 
the  site  and  time  of  the  impact. 


The  first  program  takes  the  assumed  functional  form  of  the  data  and  determines 
the  values  of  the  adjustable  parameters  which  minimize  the  error  between  the  function 
and  the  data.  The  function  used  was: 


y  «  0  for  /  <  t0 


y  —  ( mS "  —  nSm)/(m  -  rt)  for  t  >  t0 


where 

S  **  ^max  “  t)/(rmax  ~  t0) 

This  function  is  recognizable  as  Eq.  5.2c  in  which  y  is  the  pulse  amplitude,  ymax  is  its 
maximum  value,  and  .'max  is  the  time  at  which  the  maximum  occurs.  These  values  are 
fixed  by  reference  to  the  data.  The  computer  determines  the  values  of  n ,  m,  and  the 
pulse  start  time  t0  that  best  fit  the  data.  The  procedure  requires  that  initial  estimates 
of  n,  m,  are  to  be  supplied.  However,  the  answers  obtained  in  some  cases  depended 
on  the  initial  estimates.  Furthermore,  a  range  of  combinations  of  n  and  m  led  to 
equally  good  fits,  as  judged  by  the  total  RMS  error.  For  example  by  treating  n  as  a 
fixed  quantity  the  residual  error  at  the  best  values  of  t0  and  m  could  be  determined. 
This  error  was  relatively  unchanged  if  n  was  varied  by  as  much  as  several  tenths.  A 
value  of  n  —  4  for  the  data  recorded  at  the  15  cm  gage  and  n  —  6  3  for  the  data  from 
the  other  location  appeared  to  be  near  the  center  of  gravity  of  n  values  that  satisfied  to 
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error  minimization  criterion.  These  values  of  n  are  compatible  with  those  deduced  by 
the  graphical  scheme  based  on  Eq.  2.5,  i.e.  finding  an  inverse  power  of  strain  that 
behaves  linearly  with  time  at  small  strain  levels.  Furthermore,  these  values  of  n  are 
also  compatible  with  n  =  1.5  at  the  site  of  the  impact.  While  this  latter  condition  is 
not  an  absolute  criterion  it  was  the  behavior  noted  in  the  two-bar  collision  experi¬ 
ments.  In  any  case,  variations  in  the  exponent  of  a  few  tenths  in  the  observed  signal 
are  not  experimentally  significant. 

Therefore,  the  procedure  finally  adopted  was  to  fix  the  values  of  n  at  4  and  6.3  for 
the  two  gage  positions,  and  to  determine  the  value  of  m  and  /„  by  the  error  minimiza¬ 
tion  procedure. 

The  second  computer  program  uses  these  values  of  /max,  ^max,  t0 ,  «,  and  m  along 
with  the  time  t ^  of  the  first  zero  crossing  and  incorporates  them  into  Eqs.  5.12b, 
5.13,  5.14,  and  then  adds  back  in  the  residual  error  terms.  This  procedure  results  in 
the  deduced  shape  of  the  initial  pulse  in  terms  of  a  peak  normalized  strain.  The  abso¬ 
lute  strain  values  are  separately  determined  by  application  of  the  conservation  of 
momentum  requirement  Eq.  5.1  as  discussed. 

S.5  Results  of  Pulse  Correction  Procedure 

The  procedures  described  in  5.4  were  applied  to  the  pendulum  impact  data  as  well 
as  to  the  gas  gun  impact  data.  There  were  two  independent  techniques  for  estimating 
the  arrival  time  i0  of  the  pulse  leading  edge.  The  one  based  on  the  nth  root  plotting 
procedures  discussed  in  5.1  is  designated  as  t0ig  and  is  determined  by  graphical  extra¬ 
polation.  The  other  relies  on  the  best  least  square  error  fit  to  Eq.  5.2c  which  has  also 
been  described  and  leads  to  the  estimate  to  c.  Each  approach  has  limitations.  Table  5B 
summarizes  the  parameters  that  characterize  the  pulses  as  observed  prior  to  correction. 
There  is  no  pendulum  data  for  the  30  cm  gage  as  that  gage  had  not  yet  been  installed 
at  the  time  of  those  measurements.  The  pulse  maximum  relative  to  both  tog  and  tos 
are  included. 

In  order  to  test  how  well  the  correction  procedure  was  performing,  it  was  first  ap¬ 
plied  to  the  data  from  the  30  cm  gage  and  used  to  deduce  the  expected  pulse  shape  at 
the  15  cm  gage  position.  This  deduced  translated  pulse  could  then  be  compared  with 
that  experimentally  detected  by  the  15  cm  gage.  These  results  are  given  in  Table  5C. 
Two  representative  plots  which  compare  the  calculated  signal  with  the  observed  one 
are  given  in  Figs.  5.5  A  and  B.  These  curves  are  somewhat  misleading  in  that  the  pro¬ 
cedure  actually  matches  the  curves  at  the  position  of  their  peaks,  whereas  in  the 
figures  the  position  of  the  leading  edges  arc  made  to  coincide.  However,  this  way  of 
plotting  emphasizes  how  well  the  pulse  periods  match.  It  is  the  period  that  is  especial¬ 
ly  important  to  the  analysis  of  the  data. 

Finally,  the  pulse  data  from  the  15  cm  gage  position  is  translated  to  the  site  of  the 
impact  by  the  computational  procedure.  These  results  are  given  in  Table  5D.  The 
figures  of  the  corrected  pulses  are  presented  in  Chapter  6  in  conjunction  with  the  com¬ 
putations  based  on  the  corrected  pulse,  Meaningful  discussion  of  the  deduced  pulse 
shapes  is  not  possible  until  the  derived  responses  as  a  function  of  penetration  depth 
have  been  determined. 
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Table  SB 

PARAMETERS  CHARACTERIZING  OBSERVED  PULSE 
AT  15  cm  AND  30  cm  GAGES 


Specimen 

Designation 

Specimen 

Material 

Impact 

Velocity 

15  cm 

sue 

30  cm  gage 

<'•  "  Vc> 

^00  ~  'a,!* 

n 

B 

<'•  -  'o.,' 

n 

m 

180PI 

PMMA 

126.2 

mrm. 

■51 

47.98 

4 

4.85 

69.8 

6.35 

190P1 

PMMA 

I2S.I 

47.95 

4 

71.0 

111 

6.50 

190P2 

PMMA 

219.9 

34.2 

32.26 

48.5 

4 

5.11 

69.4 

6.3 

6.36 

190P3 

PMMA 

21S.7 

33.7 

31.17 

47.75 

4 

4.39 

70.8 

6.3 

6.30 

190C1 

PC 

124.5 

41.1 

40.64 

55.66 

4 

59.8 

58.1 

6.3 

6.81 

190C2 

PC 

125.9 

40.9 

40.19 

55.77 

4 

7.82 

59.0 

56.6 

i 

6.3 

6.35 

190C3 

PC 

221.5 

42.6 

39.78 

57.87 

4 

6.52 

65.5 

58.0 

6.3 

6.48 

190C4 

PC 

219.9 

41.95 

40.18 

56.96 

4 

7.11 

61.3 

57.1 

jj  ■ 

6.3 

6.30 

827P1 

PMMA 

2.0S 

153 

4 

4 

S27P2 

PMMA 

2.08 

154.5 

4 

4 

1 

827CI 

PC 

2.08 

mi 

4 

4 

jjj'.;  1 

827C2 

PC 

2.08 

wm 

4 

4 

mmm 

f*  -  time  it  pulse  maximum 
io  t  _  pulse  arrivil  time  determined  graphically 
l0(  -  pulse  arrival  time  uaing  computer  error  minimization  procedure. 


Table  5C 

COMPARISON  OF  OBSERVED  PULSE  PARAMETERS  AT  15  cm  GAGE 
WITH  VALUES  COMPUTED  FROM  DATA  AT  30  cm  GAGE 


Pulse  Parameters 


Specimen 

Designation 

^0,C 

^O.C 

m 

Calc. 

Obs. 

.  C»lc- 

Obs. 

Calc. 

Obs. 

180P1 

34.1 

32.2 

46.9 

48.0 

4.16 

4.85 

190P1 

34.6 

32  4 

47.5 

48.0 

4.29 

5.05 

190P2 

34,3 

32.3 

47.3 

48.5 

4.18 

5.11 

190P3 

34.1 

31.2 

47.3 

47.8 

4.12 

5.39 

190C1 

39.3 

40.6 

52.7 

55.7 

4.90 

7.90 

190C2 

37.8 

40.2 

50.9 

55.8 

4.46 

7.82 

190C3 

39.2 

39.8 

53.5 

57.9 

4.66 

6.52 

190C4 

38.3 

msM 

52.2 

57.0 

4.46 

7.11 
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Table  SD 


DEDUCED  PULSE  PARAMETERS  AT  SITE  OF  IMPACT 
CALCULATED  FROM  DATA  AT  15  cm  GAGE 


Specimen 

Designation 

Impact 

Velocity 

Impact  Parameters 

f-tc 

n 

m 

180P1 

126.2 

10.95 

20.95 

1.5 

1.48 

190P1 

128.1 

11.13 

21.15 

1.5 

1.58 

190P2 

219.9 

11.04 

21.41 

1.5 

1.58 

190P3 

218.7 

9.95 

20.37 

1.5 

1.18 

190C1 

124.5 

19.42 

30.20 

1.5 

4.76 

190C2 

125.9 

18.97 

29.76 

1.5 

4.60 

190C3 

221.5 

18.56 

30.05 

1.5 

3.74 

190C4 

219.9 

18.96 

30.29 

1.5 

4.18 

827P1 

2.08 

130.2 

• 

1.5 

7.7 

827P2 

2.08 

131.7 

- 

1.5 

7.7 

827C1 

2.08 

201.4 

_ 

1.5 

8.9 

827C2 

2.08 

20^.4 

" 

l.S 

8.9 
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6.  DERIVED  RESPONSE  OF  PMMA  AND  PC  TO  FREE  FLIGHT  IMPACTS 
The  corrected  pulses  can  be  kinematically  analyzed,  using  the  procedures  given  in 
the  Appendix  to  deduce  the  penetration  distance  of  the  projectile  into  the  target  as  a 
function  of  time.  Once  this  relationship  is  known,  computation  of  the  resisting  force 
as  function  of  depth,  time,  or  instantaneous  velocity  among  other  qualities  can  be  per¬ 
formed  easily.  The  analysis  also  allows  the  energy  contained  in  the  propagating  pulse 
to  be  independently  calculated.  Thus,  by  difference  the  energy  that  remains  local  to 
the  impact  can  be  calculated.  These  results  can  then  be  compared  with  the  prior 
mechanical  test  machine  data,  as  well  as  with  theoretical  and  intuitive  expectation. 
The  derived  results  for  PMMA  and  PC  follow. 

6.1  Analysis  of  PMMA  Data 

The  analysis  of  the  results  of  the  pendulum  impacts  at  2.06  m/s  and  of  the  gas  gun 
impacts  at  125  m/s  and  220  m/s  are  given  in  Figs.  6.1  A-F.  In  order  to  show  the  de¬ 
gree  of  reproducibility  the  results  of  the  duplicate  runs  are  superimposed.  Some  of 
the  major  characteristics  are  summarized  in  Table  6A.  The  circles  in  Fig.  6.1  B  are  the 
values  observed  at  5  mil  (.13  mm)  penetration  depth  in  our  prior  work  using  the 
electro-hydraulic  test  machine  at  a  constant  indentation  velocity  of  .25  and  of  2.5  m/s. 
In  Fig.  6.1  D  and  F  the  curves  obtained  in  this  prior  work  are  shown.  In  addition  the 
force-penetration  behavior  is  shown  in  Fig.  6.1  G,  in  which  the  results  for  the  three 
velocities  are  compared.  The  prior  result  obtained  at  the  2.5  m/s  using  the  electro- 
hydraulic  machine  are  included  for  reference. 

Table  6A 

SUMMARY  OF  IMPACT  RESULTS  FOR  PMMA 


Run  No. 

Velocity 

K 

£ 

Max 

Depth 

Max  Force 

%  Local  Work 

m/s 

105  lb/in. 1/2 

GPa/m1 

mils 

mm 

lbs. 

kg 

Recovered 

827P1 

2.08 

3.65 

5.24 

6.7 

.170 

155 

70.5 

56.1 

827P2 

2.08 

3.52 

5.05 

6.8 

.173 

155 

70.5 

56.0 

2.08 

3.59 

5.15 

6.8 

.172 

155 

70.5 

56.1 

180P1 

126.1 

4.30 

6  17 

37.9 

.963 

1220 

555 

26.0 

190P1 

128.1 

4.19 

6.01 

39.4 

1.000 

1210 

550 

24.7 

127.1 

4.25 

6.09 

38.7 

0.982 

1215 

553 

25.4 

190P2 

219.9 

3.28 

4.71 

71.1 

1.803 

1900 

864 

14.7 

190P3 

218.7 

5.57 

8.00 

64.3 

1.633 

1910 

868 

14.1 

219.3 

4.43 

6.36 

67,7 

1.718 

1905 

866 

14.4 

The  data  and  results  for  the  duplicate  pendulum  impact  runs  are  almost  identical. 
The  greatest  discrepancy  between  duplicate  runs  is  for  the  220  m/s  experiments.  The 
reason  for  the  differences  in  the  latter  case  is  traceable  to  the  1.09  ns  difference  in 
their  measured  periods  (the  time  interval  between  the  pulse  leading  edge  and  the  max¬ 
imum).  Such  a  difference  in  the  measured  period  is  to  be  expected  when  one  consid¬ 
ers  that  the  data  is  sampled  at  0.5  /xs  intervals  and  that  the  start  of  the  observed  pulse 
is  buried  in  the  noise  signal.  This  example  illustrates  the  importance  of  being  able  to 


define  the  start  of  the  pulse  signal.  Because  the  period  is  much  longer  in  the  case  of 
the  pendulum  impacts,  an  error  of  this  magnitude  has  a  smaller  relative  effect. 

The  force-displacement  curves  in  Figs.  6. 1  D  and  F  are  bracketed  by  the  prior  test 
machine  results  at  0.2S  and  2.S  m/s,  and  the  present  results  at  2.08  m/s  in  Fig.  6.1  B 
lie  below  the  prior  results  at  the  S  mil  penetration  check  point.  The  possible 
significance  is  discussed  in  6.3.  Fig.  6.1  G  gives  little  indication  of  a  velocity  depen¬ 
dence,  within  experimental  error.  This  is  confirmed  in  Table  6A,  in  which  the  propor¬ 
tionality  constant  K0  of  Eq.  1.1,  and  the  derived  Young's  modulus  E  are  displayed.  A 
possible  slight  trend  towards  increasing  stiffness  with  increasing  velocity  is  calculated. 
However,  there  is  a  marked  decrease  with  increasing  impact  velocity  in  the  amount  of 
energy  returned  to  the  projectile  on  rebound  from  energy  that  was  imparted  to  speci¬ 
men  in  the  locale  of  the  impact. 

6.2  Analysis  of  the  PC  Data 

The  results  for  the  impact  measurements  on  PC  are  treated  in  parallel  to  those  of 
PMMA  and  are  shown  in  Figs.  6.2  A-G  and  in  Table  6B.  The  general  features  of  the 
two  materials  are  similar.  However,  in  the  case  of  PC  there  appears  to  be  a  progres¬ 
sive  decrease  in  stiffness  as  the  velocity  increases.  Figure  6.2  B  shows  that  the  force  at 
the  S  mil  depth  value  in  the  pendulum  impact  measurements  at  2.08  m/s  is  close  to 
that  observed  using  the  test  machine  at  0.25  m/s.  Again  there  is  a  greater  discrepancy 
relative  to  the  machine  result  at  2.S  m/s.  However,  the  results  for  the  12S  m/s  and 
the  220  m/s  impacts  show  increasing  departure  from  the  machine  base  line  results. 


Table  6B 

SUMMARY  OF  IMPACT  RESULTS  FOR  PC 


Run  No. 

Velocity 

m/s 

105  lb/in.3/3 

E 

GPa/m3 

Max  Depth 

Max  Force  j 

%  Local  Work 
Recovered 

i  mils 

mm 

lbs. 

kg 

827C1 

2.08 

1.46 

2.10 

10.3 

.262 

114 

52.5 

53.0 

827C2 

2.08 

1.27 

1.82 

10.9 

.277 

111 

50.5 

52.1 

2.08 

1.37 

1.96 

10.6 

.270 

113 

51.2 

52.6 

190C1 

124.5 

0.42 

0.61 

70.6 

1.793 

877 

399 

19.8 

190C2 

125.9 

0.44 

0.63 

■ 

70.1 

1.781 

890 

405 

18.7 

125.2 

0.43 

0.62 

70.4 

1.787 

884 

402 

19.3 

190C3 

221.5 

0.25 

0.36 

126 

3.20 

1390 

632 

10.3 

190C4 

219.9 

0.31 

0.45 

128 

3.25 

1400 

636 

10.9 

220.7 

0.28 

0.41 

127 

3.23 

1395 

634 

10.6 

This  is  especially  evident  in  Fig.  6.2G. 

The  table  shows  the  same  trend  in  the  computed  values  of  K0  and  £  In  fact  the 
calculated  Young’s  modulus  at  220  m/s  is  almost  a  factor  of  5  less  than  that  of  2  m/s. 
This  decrease  of  stiffness  is  reflected  in  the  large  calculated  maximum  penetration 
depth.  Similar  to  PMMA  there  is  a  large  decrease  in  the  amount  cf  the  local  work 
recovered  as  the  impact  velocity  is  increased.  Finally  the  table  also  shows  the  general¬ 
ly  good  agreement  between  duplicate  runs. 
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6.3  Discussion  of  Results 

The  results  show  that  neither  PMMA  nor  PC  show  any  particular  trends  towards 
increased  brittleness  as  should  be  evidenced  by  an  increase  in  modulus  and  a  decrease 
in  the  viscous/ anelastic  contribution.  The  PMMA  specimens  showed  radial  crack  for¬ 
mation  at  220  m/s,  but  not  at  125  m/s.  There  is  nothing  evident  in  the  curves  to  sug¬ 
gest  that  such  a  threshold  should  occur  at  that  location.  In  fact  there  is  a  clear  trend 
with  increasing  velocity  in  both  polymers  to  retain  increasingly  more  of  the  local  ener¬ 
gy,  rather  than  to  release  it  elastically  to  the  projectile.  Hence,  instead  of  non-elastic 
processes  being  frozen  out  as  the  time  scale  decreases,  these  processes  appear  to  be 
enhanced. 

The  behavior  of  PC  in  actually  displaying  a  decrease  of  stiffness  with  increasing  im¬ 
pact  velocity  is  even  more  unexpected.  One  possibility  is  that  local  heating  due  to  the 
impact  may  soften  the  material.  However,  this  can  be  dismissed  as  unlikely  by  the  fol¬ 
lowing  argument: 

The  projectile's  kinetic  energy  of  which  90%  remains  with  the  specimen  is  9J 
at  220  m/s.  The  penetration  is  about  3.2  mm  so  that  a  crater  volume  of 
0.04  cm3  is  produced  and  a  plastic  zone  of  10  times  this  volume  results. 
Thus,  the  maximum  temperature  rise  over  this  volume  is  about  20  °  C  if  all 
of  the  energy  were  converted  to  heat.  Even  this  extreme  estimate  is  not 
sufficient  to  heat  the  PC  to  its  glass  transition.  Thus,  only  secondary  effects 
due  to  a  much  smaller  thermal  rise  are  likely  to  occur. 

Another  possibility  is  that  the  non-linear  stress-relaxation  process  (see  Eq.  1.5) 
may  dominate  over  the  usual  expected  kinetic  dependence.  The  stress-relaxation 
model  is  considered  in  the  next  chapter.  In  any  case  whereas  PMMA  is  displaying 
“normal”  behavior  based  on  conventional  wisdom,  PC  seems  to  be  behaving 
anomalously,  but  in  a  mode  favorable  to  good  impact  resistance. 

Finally,  we  consider  the  possible  reasons  why  our  present  results  give  forces  at  a 
given  depth  that  are  less  than  those  previously  measured  at  2.7  m/s  on  the  electro- 
hydraulic  equipment.  This  comparison  may  indicate  a  real  systematic  difference  be¬ 
tween  the  two  methods  such  as  due  to  the  methods  of  supporting  the  specimens. 
Another  possibility  is  that  “ringing”  or  some  other  phenomenon  in  the  electro- 
hydraulic  sensing  systems  may  be  responsible.  At  2.5  rms  that  equipment  was  being 
used  at  the  upper  limit  of  its  capability,  where  it  is  known  that  such  complications  can 
begin  to  have  an  effect.  The  data  that  we  had  obtained  on  that  machine  showed  a 
somewhat  larger  than  expected  increase  in  stiffness  when  the  velocity  was  increased 
from  0.25  m/s  to  2.5  m/s  relative  to  the  data  at  0.25  m/s  and  less.  To  what  extent,  if 
any  the  source  of  the  discrepancy  is  attributable  to  the  electro-mechanical  machine 
would  have  to  be  resolved  by  future  experimentation.  In  any  case  the  present  results 
are  in  good  agreement  with  the  prior  one  at  0.25  m/s.  It  is  also  important  to  keep  in 
mind  that  velocity  is  not  constant  in  the  present  experiments  so  that  the  results  are 
not  directly  comparable.  However  as  Fig.  6.3A  shows  for  the  example  of  the  220  m/s 
impact  velocity  using  a  PC  target,  the  velocity  remains  constant  to  within  10%  over  the 
first  half  of  the  penetration. 
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7.  VISCO-ELASTIC  MODEL  OF  POLYMER  IMPACT  BEHAVIOR 

The  m^jor  features  of  the  polymer  response  noted  in  this  study  indicates  (1)  the 
maximum  reaction  force  precedes  the  maximum  penetration  depth,  (2)  the  change  of 
load  with  penetration  distance  is  more  pronounced  during  projectile  exit  than  on  entry, 
(3)  the  load  on  exit  drops  to  zero  when  the  penetration  depth  is  an  appreciable  frac¬ 
tion  of  the  maximum  depth,  and  (4)  most  of  the  local  energy  absorbed  by  the  impact 
is  retained  at  the  time  of  exit.  To  this  behavior  should  be  added  such  prior  observa¬ 
tions  as  (5)  yielding  densities t ion  occurs  during  the  projectile  entry  under  the  impact 
site,  (6)  a  strain-rate  and  history-dependent  stress  relaxation  accompanies  the  strain 
process,  and  (7)  over  many  decades  of  velocity  there  is  a  trend  towards  a  net  harden¬ 
ing  with  increased  penetration  rate.  The  latter  observation  appears  to  have  been  re¬ 
versed  by  PC  at  penetration  rates  of  about  2  m/s  and  greater. 

Even  though  the  separate  phenomena  themselves  are  complex,  it  is  instructive  to 
consider  how  they  may  be  interrelated.  Accordingly  a  relatively  simple  conceptual 
model  has  been  formulated  to  try  to  connect  the  processes  (5) -(7)  so  as  to  qualitative¬ 
ly  predict  or  explain  the  present  results  (l)-(4).  The  model  uses  simple  mechanisms 
with  the  objective  of  gaining  insight  into  the  individual  differences  in  the  impact 
behavior  of  PC  and  PMMA. 

Briefly  the  penetration  is  assumed  to  result  in  elastic  and  viscous  deformation.  The 
viscous  flow  that  occurs  is  ultimately  recoverable  upon  annealing  and,  therefore,  is 
really  an  anelastic  response.  The  polymers  are  conceived  as  being  a  collection  of  in¬ 
tertwined,  snarled  springs  embedded  in  a  viscous  matrix.  The  viscous  and  the  elastic 
reaction  forces  to  an  applied  load  are  partly  in  parallel  and  partly  in  series.  As  defor¬ 
mation  occurs,  accompanied  by  some  alignment,  the  mix  of  parallel  to  series  response 
can  change.  Finally,  the  viscosity  of  the  matrix  itself  is  history  dependent,  but  is  oth¬ 
erwise  Newtonian.  The  resistance  to  the  sphere  penetration  is  partly  due  to  a  Stokes- 
like  viscous  resistance,  and  partly  to  a  Hertz  elastic  response.  This  leads  to  a  net  force 
vs.  penetration  relationship  which  in  turn  is  velocity  and  time  dependent.  The  details 
and  results  are  discussed  next. 

7.1  Details  of  the  Model 

The  elastic  and  viscous  response  is  considered  to  be  connected  by  the  spring- 
dashpot  system  shown  below.  This  system  differs  from  the  more  conventional 
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models  in  that  the  left  hand  series  arrangement  and  the  right  hand  parallel  arrange¬ 
ment  are  connected  by  a  hinged  bar  so  that  each  sub-system  contributes  separately  to 
the  total  response.  By  placing  the  compressive  load  a  fractional  distance  /along  the 
bar,  the  response  will  be  /  parts  parallel  and  (1-./)  parts  series.  The  two  dashpots 
represent  the  viscous  response  of  the  same  polymeric  material,  and  are,  therefore, 
assumed  to  behave  identically.  Similarly  the  springs  are  assumed  really  to  be  a  single 
spring  which  in  total  force  is  Hertzian  with  respect  to  total  elastic  penetration,  but  the 
force  is  split  between  the  two  modes  of  material  behavior. 

In  following  expressions,  the  subscripts  refer  to  the  system  element  defined  by  the 
figure.  The  subscript  zero  refers  to  the  quantities  as  measured  at  the  point  of  applica¬ 
tion  of  the  force,  i.e.  at  the  projectile.  Force  is  F,  displacement  x,  and  the  time  rate 
of  displacement  is  given  by  the  superscripted  dot  notation.  The  variables  are  F,  x, 
and  x  at  each  of  the  five  locations,  so  that  IS  equations  are  needed  to  define  the 
values.  However,  by  definition  F,  -  F2,  x3  -  x4,  and  x3  -  x4,  and  we  can  consider 
x0  and  x0  as  known,  from  the  solution  of  the  equation  of  motion  for  the  projectile. 
Therefore,  10  equations  suffice.  Designating  the  dashpot  resistance  by  n,  the  spring 
constant  by  k  and  the  fraction  parallel  by  /,  these  are: 


(a) 

F  -  F,  +  F3  +  F4 

( b ) 

fF\  -  (1  -  /)  (F3  +  F4) 

(c) 

x  -  (1  -  f)  (x,  +  x2)  +  /x3 

id) 

X  -  (1  —  /)  (xt  +  x2)  +  /x3 

(e) 

F7  - 

(/) 

F4  -  fix 4 

(8,h,i) 

X,.  -  J*  x,dt  i  -  1,2,4 

CO 

F,  +  F3  -  k[(1  -  /)x,  +  /x313/2 

(*) 

F,  -  kI(1-/)x,]3/2 

or 

=  (y)it(l  -  f)xl/2X  1  for  /=  1 

Equation  (a)  is  simply  the  balance  of  forces;  (ft)  is  the  balance  of  moments  and  is 
needed  to  satisfy  conservation  of  energy;  (c)  and  (d)  are  geometric  requirements;  (<?) 
and  (J)  are  the  dashpot  responses;  (g,h,i)  are  obvious;  (/)  represents  the  force  of  the 
non-linear  Hertzian  spring  in  which  the  total  elastic  force  is  assumed  to  be  determined 
by  the  total  elastic  displacement  as  measured  at  the  ball;  (k)  is  an  approximate  expres¬ 
sion  for  the  series  elastic  force.  The  constant  *  can  be  identified  with  K0  of  Eq.  1.1. 


The  dashpot  constant  n  can  be  defined  via  the  hydrodynamic  law  for  the  motion  of 
a  sphere  through  a  viscous  medium: 


Force  «  3  UxA/lRQ 


(Eq.  7.2) 
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in  which  ^  is  fluidity,  U  the  velocity  of  the  sphere  of  radius  R  relative  to  the  fluid  and 
A  the  surface  area,  and  the  force  is  applied  in  the  direction  of  motion.  This  equation 
is  valid  for  portions  of  a  sphere  assumed  to  be  in  an  infinite  medium.  For  the  present 
case,  the  term  R  can  be  interpreted  as  a  boundary  layer  and  will  be  replaced  by  the 
distance  from  the  ball  surface  to  the  boundary  of  the  plastic  or  densified  zone.  Since 
H  is  force/ velocity,  and  since  A  -  2irRx,  Eq.  7.2  can  be  written  as: 

H  —  3irRx/QRT  (Eq.  7.3) 

where  Rt  is  the  boundary  layer  thickness.  The  fluidity  is  strain-rate  dependent  and 
will  be  treated  later. 

The  computational  procedure  is  to  assign  the  proper  material  parameters  to  k  and 
fi ,  assume  a  value  of  /,  given  an  initial  impact  velocity.  At  each  time  increment  an 
updated  value  for  x,  and  /x  are  computed.  The  above  system  of  equations  are  solved 
to  finally  yield  the  reaction  force  at  that  instant  of  time.  This  determines  the  decelera¬ 
tion  and  allows  the  next  increment  of  x  to  be  computed.  As  this  process  continues, 
the  projectile  stops  and  changes  direction,  being  propelled  by  the  residual  energy  in  the 
springs.  During  this  process,  the  dashpot  in  the  series  arm  is  constantly  relaxing,  so 
that  this  portion  of  the  energy  is  not  recovered.  Similarly  the  dashpot  in  the  parallel 
arm  inhibits  the  release  of  energy  from  its  spring  so  that  not  all  of  that  stored  energy 
is  available  to  the  projectile  during  its  rebound  trajectory.  It  should  be  noted  that  this 
computational  scheme  is  approximate  and  that  the  results  can  exhibit  artificial  irregu¬ 
larity. 

7.2  The  Time  Dependence  of  the  Fluidity 

In  our  prior  work  a  model  was  proposed  to  account  for  the  following  experimental 
load  relaxation  law 

Fit)  -  F(0)(1  +  At)~B  (Eq.  1.5) 

This  law  was  observed  when  the  penetration  rate  was  abruptly  changed  from  a  cor- 
stant  velocity  to  zero.  In  addition,  the  model  qualitatively  accounted  for  the  increases 
and  decreases  in  compliance  when  the  penetration  rates  were  abruptly  changed. 

The  elements  of  the  model  are: 

1.  There  is  a  residual  intrinsic  fluidity  4>0  in  the  fully  annealed,  unstrained  polymer. 

2.  Fluidity  can  be  increased  by  increasing  the  volume  concentration  of  unpinned 
sites.  The  polymer  is  envisioned  as  being  tangled,  and  that  unsnarling,  or  freeing 
the  polymer  chain  allows  relative  motion  to  occur. 

3.  The  unpinning  is  facilitated  by  straining  the  polymer,  which  is  expected  to  align 
the  molecular  segments  by  stretching  or  shearing. 

4.  Unpinning  is  a  cooperative  process  and  is  made  easier  by  prior  unpinning. 

5.  There  is  a  spontaneous  tendency  of  the  molecules  to  rekink.  This  is  a  first  order 
reaction  and  the  rate  is  proportional  to  the  mobility  of  the  system,  i.e.  to  fluidity. 

This  results  in  the  following  system  of  equations: 

<J>  -  d>0  +  (Eq.  7.4) 

i  -  A«((l  -  o  -  e<D{  (Eq.  7.5) 

in  which  n  is  the  volume  density  of  potentially  unpinnable  sites,  0  is  the  fluidity  con- 
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tribution  per  site,  £  is  the  fraction  of  unpinned  sites,  b  is  a  constant  reflecting  the 
efficiency  of  the  unpinning  process,  c  is  a  rate  constant  for  the  decay  of  the  pinned  site 
population,  and  e  is  the  applied  strain  rate. 

Below  the  glass  transition  temperature,  one  expects  to  be  very  small.  There¬ 
fore,  taking  the  time  derivative  of  Eq.  7.4  and  inserting  7.5  gives 

1  -  £)  -  cn<f»i2  (EQ*  76a> 

-  a4»  -  /3$2  (Eq.  7.6b) 

in  which 

a  tm  bf  (Eq.  7.6c) 

and 

/3  -»  c  +  (be/ n<f>)  (Eq.  7.6d) 

Finally  we  can  relate  i  to  x2  and  x3  in  the  spring-dashpot  model  of  the  preceding 
section,  by  «,  -  x,/R2 . 

At  steady  state  -  0,  or 

-  a/fi  (E{1-  7  7> 

The  values  for  the  constants  can  be  determined  from  the  experimental  values  ob¬ 
served  for  the  load  relaxation  law. 

For  the  case  of  a  pure  series  arrangement,  it  was  previously  shown  that 
F-  Ke 

- 

in  which  K  —  AE  and  J  —  3A/2.  This  finally  led  to 
c  -  2E/3B 

in  which  Band  care  the  constants  in  Eq.  1.7,  and  7.6d  respectively,  and  E is  Young’s 
modulus.  It  can  be  shown  that  as  long  as  the  fraction  of  unpinned  sites  to  pinned  sites 
is  small,  as  is  expected  to  be  the  case,  then  the  term  in  parentheses  in  Eq.  7.6d  can  be 
neglected,  or 

f}  ~  2E/3B  (Eq.  7.8) 

Similarly  the  prior  work  showed  that  under  the  same  approximate  conditions, 

A  -  d>„c 

which  upon  introduction  of  7.8  and  7.7  leads  to 

a  ~  A,  (Eq.  7.9) 

in  which  A  is  an  experimentally  measured  quantity,  itself  a  function  of  indentation  ve¬ 
locity. 

These  relationships  can  then  be  introduced  into  the  computation  scheme  outlined 
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in  the  previous  section.  In  this  way,  by  introducing  the  dependence  of  A  on  velocity, 
and  through  the  use  of  Eq.  7.6b,  the  fluidity  is  made  history  and  time  dependent. 

7.3  Discussion  of  Input  Data  and  Procedure 

In  order  to  carry  out  the  computations  specific  values  must  be  designated  for  the 
Young’s  modulus  E,  the  ratio  X  of  the  densified  zone  relative  to  the  displaced 
volume,  the  values  of  A  and  B  appropriate  for  the  relaxation  law  Eq.  1.7,  and  the  frac¬ 
tion  /of  parallel-like  behavior  in  Eq.  7.1a-j. 

Previously,  the  size  of  the  densified  zone  in  PC  was  found  to  be  independent  of 
velocity.  We  assume  this  is  also  the  case  for  PMMA.  Therefore,  those  values  of  X 
for  PC  and  PMMA  of  3.6  and  6.2  respectively  were  introduced  into  the  model  for  the 
220  m/s  impact  by  a  4.5  mm  steel  ball. 

In  our  prior  work  the  observed  relaxation  behavior  was  studied  in  detail  using  the 
mechanical  test  machines.  These  results  showed  that  B  is  relatively  insensitive  to  ve¬ 
locity  u,  whereas  A  is  nearly  proportional  to  u  Therefore,  it  is  convenient  to  define 
A0  as 

A0  -  A/u 

where  A0  has  the  dimensions  of  an  inverse  length,  but  may  be  a  weak  function  of  ve¬ 
locity.  Both  A0  and  B  are  functions  of  depth  as  well.  The  values  for  A0  and  B  were 
obtained  for  PMMA  over  a  5  decade  span  of  velocity.  However,  that  range  was, 
nevertheless,  4  decades  of  velocity  below  that  of  our  current  work.  The  corresponding 
data  for  PC  is  less  complete  and  exhibits  more  scatter.  However,  there  appears  to  be  a 
regular  trend  in  the  product  (AB)  to  increase  approximately  linearly  with  the  loga¬ 
rithm  of  velocity.  Because  it  is  the  product  AB  that  defines  the  steady  state  fluidity, 
we  somewhat  arbitrarily  fixed  B  at  the  experimental  values  obtained  at  low  velocities 
and  considered  A0  to  be  velocity  dependent.  The  values  of  B  used  in  the  computa¬ 
tions  were 

B  -  0.05  +  0.05yfX[R  for  PMMA 


and 

5-  0.025  +  0.03 (X/R)  for  PC 

where  the  approximate  dependence  on  penetration  depth  X  is  taken  into  account. 

The  Young’s  modulus  E  is  also  expected  to  be  velocity  dependent.  One  expects  a 
priori  that  E  should  increase  slowly  with  increasing  velocity  to  some  ultimate  upper 
bound  plateau.  Apparent  values  are  available  from  the  results  of  the  previous  chapter. 
In  the  case  of  PC  the  velocity  dependence  of  the  calculated  modulus  was  anomalous. 
Therefore,  it  was  of  interest  to  find  out  if  the  behavior  of  PC  could  be  attributable  to 
the  velocity  dependence  of  A0  or  to  the  balance  between  series  and  parallel  behavior 
as  measured  by  f. 

The  approach  was  to  consider  X  and  B  as  fixed  known  quantities,  and  A0,  E,  and  / 
as  quantities  to  be  determined  from  the  experimental  pulse  data  for  PC  and  PMMA. 

Ideally,  the  pulse  computed  by  this  model  should  match  the  experimental  pulse 
everywhere.  In  practice,  the  most  important  details  to  be  matched  are  the  initial  force 
dependence  at  the  start  of  the  pulse,  the  magnitude  and  time  position  of  the  peak,  and 
the  exit  time  and  exit  velocity  of  the  projectile.  If  /,  E,  and  A0  remain  constant  dur¬ 
ing  an  impact,  the  initial  force  behavior  is  determined  mainly  by  E  Therefore,  there 
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is  little  latitude  in  adjusting  the  value  of  the  modulus.  However,  / and  An  have  largely 
compensating  effects  on  the  pulse  shape.  Increasing  /tends  to  increase  the  stiffness  of 
the  system,  to  increase  the  peak  force,  to  shorten  the  time  to  reach  the  maximum  and 
to  exit,  and  usually  to  increase  the  exit  velocity.  Similar  effects  result  from  decreasing 
A0 .  If  the  peak  force  is  held  constant  by  applying  compensating  charges  in  A0  and  /, 
then  the  effect  of  decreasing  A0  is  to  shift  the  peak  to  shorter  times,  to  reduce  the 
maximum  penetration  depth,  to  extend  the  exit  time  and  increase  the  exit  velocity, 
and  to  decrease  the  depth  of  the  impact  crater  at  exit.  This  array  of  consequences  due 
to  just  two  parameters  A0  and  /  is  not  necessarily  compatible  with  improving  the 
overall  match  between  the  modelled  result  and  the  experimental  values. 

7.4  Computational  Results 

The  expected  responses  of  PC  and  PMMA  when  impacted  by  a  4.5  mm  steel  ball  at 
220  m/s  have  been  calculated  using  the  preceding  model.  Some  representative  “best 
fits”  are  given  in  Table  7A  for  these  two  materials,  using  several  combinations  of  the 
impact  variables.  The  calculated  results  in  terms  of  penetration  depth,  and  local  ener¬ 
gy  were  plotted  in  the  same  manner  as  were  the  experimental  pulse  data  in  Figs.  7.4  A 
and  B  for  PC  and  Figs.  7.4  C  and  D  for  PMMA.  Comparison  with  the  experimental 
data  is  shown.  The  force-time  results  for  PMMA  can  be  seen  to  correspond  qualita¬ 
tively  with  the  experimental  results,  except  that  the  calculated  pulse  is  skewed  towards 
longer  times.  The  calculated  pulse  shape  for  PC  is  similarly  skewed  but  drops  precipi¬ 
tously  on  the  descending  limb.  These  curves  for  PC  and  PMMA  correspond  to  the 
columns  labelled  Case  2  in  Table  7A. 

The  most  difficult  characteristics  to  satisfy  were  the  time,  velocity,  and  penetration 
depth  at  exit.  As  the  table  shows  the  fit  in  the  case  of  PC  in  regard  to  exit  velocity  is 
rather  poor.  However,  the  computer  program  is  somewhat  unstable  in  the  descending 
region  when  /is  close  to  one.  Thus,  the  computed  exit  velocity  is  not  completely 
reliable. 

The  results  from  the  model  suggest  a  rather  different  behavior  between  PC  and 
PMMA.  The  former  has  a  much  smaller  Young’s  modulus  and  exhibits  a  nearly  pure 
parallel-like  behavior  with  a  relatively  fluid  damping  component.  The  PMMA  has  a 
more  than  10-fold  larger  Young’s  modulus,  responds  in  a  series-like  fashion  but  with 
a  more  viscous  damping  component.  These  differences  appeared  to  persist  over  all  of 
the  variations  of  the  parameters  examined  to  try  to  match  the  modelled  pulse  shape  to 
the  experimental  data. 

If  the  assumed  modulus  for  PC  were  increased  to  a  level  more  like  that  for 
PMMA,  the  fit  at  the  initial  stages  of  penetration  was  quite  poor.  However,  if  /  were 
to  be  allowed  to  decrease  with  increasing  time  or  penetration  depth,  it  may  be  possible 
to  reproduce  the  experimental  curve.  This  complicates  the  model.  However,  there  is 
no  a  priori  reason  to  accept  or  reject  this  additional  degree  of  freedom. 

The  model,  while  conceptually  simple,  is  computationally  complex.  The  computer 
solution  at  times  is  not  adequately  damped  and  takes  on  an  oscillatory  character  or 
even  instability.  This  could  be  remedied  by  developing  a  better  algorithm.  The  exit 
depth  and  velocity  were  the  most  difficult  properties  to  try  to  match,  without  severly 
compromising  the  fit  elsewhere. 

In  the  case  of  PC  a  fixed  value  of  /  does  not  reproduce  the  descending  region  of 
the  pulse  well.  However,  in  this  case  the  computer  solutions  for  high  parallel  contents 
tend  to  behave  as  noted.  Until  the  model  is  more  completely  explored  it  is  not  clear 
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whether  the  problem  in  fitting  the  exit  values  is  a  weakness  of  the  model  itself,  or  an 
indication  of  the  complexity  of  the  material  behavior.  In  the  case  of  PMMA  the  de¬ 
scending  region  is  modelled  reasonably  well  by  a  single  fixed  value  of  /  for  the  entire 
impact  process. 

The  values  of  A0  required  to  provide  the  best  matches  with  the  data  are  a  factor  of 
2-5  greater  than  would  be  expected  from  a  simple  linear  extrapolation.  (In  the  case  of 
PC  the  extrapolated  value  is  about  1500  and  for  PMMA  about  220.)  If  the  values  de¬ 
duced  from  the  model  are  significant,  then  these  should  provide  a  basis  for  a  priori 
impact  response  predictions  at  both  lesser  and  greater  velocities  than  used  in  the 
present  examples.  This  would  provide  an  addition  test  of  the  validity  and  utility  of  the 
model. 


Table  7A 

COMPARISON  OF  CALCULATED  AND  OBSERVED  PULSES  FOR  PC  AND  PMMA 
RESULTING  FROM  INPACT  WITH  4.5  mm  STEEL  BALL  AT  220  m/s 


PC 

PMMA 

Parameter 

Model 

Model 

Model 

Model 

Exptl 

Case 

Case 

Exptl 

Case 

Case 

1 

L 

1 

2 

Force  at  2  /xs 
(lbs.) 

82 

130 

104 

716 

690 

685 

Force  at  5  /*s 
(lbs.) 

324 

367 

330 

1460 

1580 

1460 

Max  Force 
(lbs.) 

1400 

1390 

1450 

2030 

2080 

1990 

Time  at 

Peak  (|is> 

17.6 

21 

20 

9.2 

10 

10 

Max  Penetration 
(Mils) 

120 

124 

121 

61.7 

63 

71 

Exit  Time 

MS 

30 

32 

* 

20 

21 

20.5 

Exit  Velocity 
m/s 

67 

121 

• 

77 

112 

86 

Depth  at  Exit 
(Mils) 

113 

94 

* 

:o 

35 

54 

A0  (cm-1) 

- 

3000 

10000 

1 

400 

3000 

E  (GP.) 

- 

.52 

.55 

6 

8.5 

/(fraction 

parallel) 

* 

.82 

.94 

0 

.25 
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8.  DISCUSSION  AND  CONCLUSIONS 

8.1  Discussion  of  Measurement  and  Analysis  Procedures 

A  way  to  treat  pulse  data  has  been  developed  to  deduce  the  details  of  the  force 
response  as  the  projectile  penetrates  the  target.  Although  the  methodology  is  derived 
from  empirical  observation,  rather  than  derived  from  the  exact  theory,  the  procedure 
is  consistent  with,  and  in  part  determined  by  basic  physical  requirements.  For  exam¬ 
ple,  the  maximum  rate  of  propagation  of  a  sonic  disturbance,  enters  into  the  formula¬ 
tion  naturally,  as  does  the  normal  rate  of  travel  of  plane  waves  in  a  long  cylinder  and 
the  way  the  leading  edge  changes  form  as  it  travels. 

The  elementary  theory,  which  assumes  a  single  characteristic  sound  velocity,  is 
adequate  and  useful  when  the  period  from  the  start  to  the  peak  of  the  pulse  is  long 
compared  with  the  spread  of  arrival  times  of  the  components  travelling  at  the  different 
velocities.  In  the  case  of  the  2.5  cm  diameter  PMMA  impact  bar  used  in  this  study, 
and  the  sensor  located  6  diameters  from  the  bar  end,  this  time  spread,  due  to  the  ve¬ 
locity  differences,  is  about  20  /as.  Hence,  if  in  a  pulse  with  a  period  of  ten  times  this 
value,  this  transient  was  not  detected,  or  was  ignored,  the  pulse  could  nevertheless  be 
analyzed  with  less  than  a  10%  error  using  the  elementary  theory.  Such  is  the  case  for 
impacts  of  velocities  below  about  2  m/s.  For  higher  velocities,  a  procedure  such  as 
that  described  in  Section  6,  is  required. 

The  experimental  measurements  are  straightforward.  Part  of  the  complications  due 
to  dispersion  could  be  mitigated  by  the  use  of  a  bar  of  some  other  material.  In  studies 
of  impact  on  materials  other  than  polymers,  this  is  a  recommended  choice.  However, 
in  the  study  of  impacts  on  polymers,  requirements  of  good  acoustic  impedance  match¬ 
ing  between  the  bar  and  the  specimen  necessitate  using  a  polymeric  bar.  While  increas¬ 
ing  the  magnitude  of  the  corrections  needed  to  deduce  the  true  originating  pulse, 
these  complications  remain  manageable. 

The  amount  of  effort  involved,  and  the  way  finally  established  to  correct  observed 
pulses  of  short  duration  was  unexpected  in  as  much  as  impact  bar  measurements  are  a 
common  experimental  procedure.  Furthermore,  the  basic  equations  governing  acous¬ 
tic  propagation  are  known  as  well  as  are  Maxwell's  equations  in  the  area  of  classical 
electricity  and  magnetism.  Nevertheless,  the  solution  of  the  problem  of  leading  edge 
behavior,  is  very  complex.  The  functional  forms  that  we  observed  experimentally  are 
recommended  to  the  theoretician  interested  in  solving  this  leading  edge  problem.  It 
may  be  that  the  pulse  shape  found  and  its  mode  of  propagation  are  peculiar  to  the  ini¬ 
tial  and  boundary  conditions  of  our  particular  experiment,  in  as  much  as  the  pulse  was 
initiated  essentially  at  a  point  on  the  center  of  the  bar  end  in  the  present  work. 

The  observed  leading  edge  transient  has  the  form 

F(x,t )  -  Kt”+a *  (Eq.  8.1) 

where  x  is  distance  from  the  site  of  impact,  K ,  n  and  a  are  constants,  and  t  is  locai 
time  measured  from  the  first  arrival  of  the  pulse.  This  transient  was  assumed  to  be 
accommodated  into  the  overall  pulse  by  the  equation 

F  —  At"  —  Btm  (Eq.  8.2) 

This  may  not  be  the  best  representation  of  the  way  the  transient  and  the  overall  pulse 
geometries  are  related.  A  better  description  would  require  further  investigation.  Our 
justification  was  the  excellent  empirical  fit  found  with  our  observed  data.  In  any  case 
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the  deviations  from  this  law  were  reinserted  into  the  final  result.  In  trying  to  deter¬ 
mine  the  evolution  of  the  pulse  shape,  a  central  assumption  and  requirement  in  our 
analysis  was  that  once  the  pulse  is  created,  its  propagation  is  only  a  function  of  the 
bar.  Thus,  any  artifacts  that  the  procedure  described  in  5.4  introduces  into  the  final 
results,  should  be  independent  of  the  target  material  or  impact  dynamics. 

8.2  Discussion  of  the  PC  and  PMMA  Results 

The  impact  results  of  the  two  materials  show  distinctly  different  trends.  One  may 
question  how  well  the  computed  results  reflect  the  true  force-time-penetration  rela¬ 
tionships  of  the  actual  impact  event.  The  results  for  PMMA  were  largely  velocity  in¬ 
dependent,  and  were  compatible  with  the  prior  results  obtained  on  the  electro- 
hydraulic  test  machine.  Those  for  PC  show  a  progressive  softening  with  increasing  ve¬ 
locity.  The  results  for  PC  and  PMMA  using  the  pendulum  impact  method  confirmed 
the  prior  test  machine  results  at  comparable  velocity.  In  the  case  of  the  gas  gun  im¬ 
pacts  the  most  sensitive  factor  affecting  the  results  is  the  value  of  the  computed  time 
period  after  shifting  the  pulse.  This  value  is  in  turn  affected  by  the  accuracy  with 
which  the  period  at  the  gage  is  determined,  and  by  the  velocity  values  used  to  effect 
the  shift  in  shape  and  position. 

The  error  in  determining  the  period  at  the  gage  is  mainly  due  to  the  error  in 
detecting  the  position  of  the  leading  edge.  The  graphical  procedure  and  the  computer 
error  minimization  procedure  often  led  to  differences  of  several  /lis.  The  error  minim¬ 
ization  was  in  principle  the  more  objective  measure,  and  was  used  to  refine  the  graphi¬ 
cal  estimate.  In  spite  of  the  high  precision  of  the  computed  result,  its  reliability  is  es¬ 
timated  to  be  ±  1  fis.  Such  an  error  is  insufficient  to  account  for  anomalous  softening 
deduced  from  the  PC  results. 

The  other  factor  affecting  the  calculated  shifted  pulse  period  is  the  difference  in  the 
sound  propagation  velocities.  Even  if  incorrect  velocity  values  were  used,  they 
remained  unchanged  in  the  analysis  process.  Therefore,  since  the  pulse  period  in 
PMMA  is  smaller  than  PC,  there  would  be  an  even  larger  artifact  produced  in  the 
former  material.  The  only  remaining  possibility,  if  the  softening  effect  were  not  real, 
is  that  the  propagation  velocity  of  the  leading  edge,  or  of  other  parts  of  the  pulse  were 
different  for  PC  than  for  PMMA.  This  cannot  be  disproved,  but  certainly  contradicts 
the  principle  that  the  pulse  propagation  in  the  bar  depends  only  on  the  bar. 

In  summary,  it  appears  that  the  increasing  impact  softening  in  PC  with  increasing 
impact  velocity  is  real.  Since  elastic  modulus  is  expected  to  increase  with  increasing 
velocity,  the  conclusion  is  that  the  viscous  response  of  PC  is  becoming  more  pro¬ 
nounced.  This  would  have  a  beneficial  effect  on  enhancing  the  impact  damage  resis¬ 
tance  of  PC  relative  to  PMMA, 

8.3  Results  of  Model  Calculations 

The  initial  force-penetration  behavior  of  PMMA,  as  well  as  the  time  and  depth  of 
the  pulse  maximum  are  well  approximated  by  assuming  the  response  to  be  predom¬ 
inantly  a  series  arrangement  of  the  elastic  and  viscous  flow  response.  However,  the 
behavior  of  PC  is  much  more  complex.  This  material  is  initially  much  softer  than 
predicted  on  the  basis  of  the  expected  continuing  tendency  of  the  elastic  modulus  to 
increase  with  increasing  velocity.  However,  with  deeper  penetration  the  response 
stiffens,  This  can  be  reconciled  with  a  predominantly  parallel  elastic-viscous  element 
grrangement  in  combination  with  a  greater  dashpot  fluidity. 
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However,  this  arrangement  results  in  a  very  steep  descending  limit  and  a  small  re¬ 
bound  velocity.  The  response  can  only  be  made  to  conform  more  with  experiment  by 
an  increasing  shift  to  an  increasingly  series-like  arrangement.  This  leads  to  the  con¬ 
cept  of  PC  behaving  like  an  elastic  skeleton  of  which  the  interstices  are  filled  with  a 
relatively  mobile  fluid.  However  as  the  strain  level  is  increased,  the  elastic  network 
breaks  up.  In  terms  of  the  model  this  requires  that  the  parameter  /decrease  with  in¬ 
creasing  strain. 

Such  a  mental  picture  would  account  for  the  decreasing  modulus  with  increasing 
strain-rate  (impact  velocity)  if  the  fraction  of  molecules  associated  with  the  elastic 
skeleton  were  to  decrease  with  increasing  strain  rate.  This  is  consistent  with  the  time 
dependent  fluidity  model,  in  which  we  considered  an  unpinning  to  control  fluidity. 
The  greater  the  strain  rate,  the  greater  the  degree  of  unpinning,  according  to  Eq.  7.4 
and  7.5.  As  long  as  the  fraction  (  of  unpinned  sites  is  much  less  than  unity,  the  un¬ 
pinning  would  have  no  effect  on  the  modulus.  However,  if  {  becomes  large,  then  the 
remaining  skeleton  is  weak  and  a  decrease  of  elastic  modulus  is  expected.  According 
to  Eqs.  7.9,  and  7.6b,  {  will  become  large  if  the  relaxation  parameter  A  -  A0x  velocity 
is  large.  In  the  case  of  PC,  the  best  estimates  for  A0  were  several  fold  greater  than 
the  values  for  PMMA. 

This  interpretation,  while  speculative,  provides  a  possible  explanation  for  the 
difference  in  behavior  of  the  two  polymers.  It  also  suggests  that  one  of  the  important 
variables  affecting  the  resistance  of  a  polymer  is  the  A0  parameter.  Thus,  it  would  be 
useful  to  understand  in  terms  of  the  molecular  structure  how  to  modify  this  property. 
Finally,  the  relaxation  model  and  the  derivative  impact  model  require  considerable 
critical  testing  before  they  can  be  considered  reliable  guides  to  the  behavior  of  real 
materials. 
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11.  APPENDIX 


11.1  Conversion  of  Impact  Pulse  Data  Into  Force-Time  and  Energy-Time 
Relationships 

The  basis  for  deducing  the  force  indentation  response  in  the  impact  event  from  the 
impact  pulse  was  discussed  in  previous  reports(5)  and  is  incorporated  here  for  conveni¬ 
ence. 

Because  the  impact  occurs  over  a  concentrated  area,  the  compressive  acoustic  pulse 
radiates  initially  as  a  spherical  wave.  However,  as  it  travels  down  the  bar,  the 
compressive  wave  becomes  increasingly  planar.  If  the  initial  impact  does  not  occur  on 
the  axis  of  symmetry  of  the  bar,  then  flexural  modes  can  be  excited  in  the  bar. 
Furthermore,  surface  and  transverse  waves  can  be  expected  as  well.  Thus,  the  signal 
detected  at  the  strain  gage  comprises  a  complex  superposition  of  waves.  In  spite  of 
this,  the  data  will  be  treated  as  if  the  strain  pulse  consisted  only  of  the  planar  compres¬ 
sive  component.  This  means  the  analysis  will  be  somewhat  in  error  and  that  calibra¬ 
tion  to  provide  the  necessary  corrections  would  be  desirable. 

For  purposes  of  analysis,  it  is  convenient  to  refer  to  the  position  of  the  struck  end 
of  the  bar  as  the  origin  of  the  z  axis  and  the  time  t  =  0  as  the  instant  of  first  contact 
between  the  ball  and  the  bar.  The  ball  has  a  mass  m  and  is  assumed  to  be  incompres¬ 
sible. 


Deceleration  of  the  spherical  projectile  results  from  the  local  opposing  material 
reaction  force  which  in  turn  is  supported  elastically  in  the  impac.  bar.  Hence,  the  in¬ 
stantaneous  deceleration  force  is  Acr(t)  where  A  is  the  cross  section  of  the  bar,  and  <r 
the  area-averaged  longitudinal  stress.  Therefore, 


na  —  —  Acr(t) 


(Eq.  a) 


(Eq.  b) 


This  can  be  integrated  directly  to  give 

z  -  -  ( AE / m)J*  edt  +  v0 
0 

and 


(Eq.  c) 

(Eq.  d) 


where  Eis  the  Young’s  modulus  of  the  bar,  epsilon  is  strain,  and  v0  is  the  initial  ve¬ 
locity  of  the  projectile.  However,  the  collision  causes  the  impact  bar  to  be  elastically 
compressed.  Therefore,  the  end  of  the  bar  is  displaced  in  the  direction  of  the  impact¬ 
ing  projectile  a  distance  w 


(Eq.  e) 
(Eq.  e') 


where  c  is  the  velocity  of  sound.  Thus,  penetration  is  given  by 

X  -  z  -  w  (Eq.  f) 

and  the  velocity  of  the  projectile  relative  to  the  end  of  the  bar  is  given  by 


S.  J;  ‘u£.L* , 


■ujTiMSWuifca': i_»  “ 
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x  -  i  -  ce 


(Eq.  g) 

Assuming  the  projectile  to  rebound  with  a  velocity  v /,  its  change  in  momentum  must 
equal  the  momentum  imparted  to  the  bar.  This  is  the  total  time  integral  of  the  force 
(neglecting  reflections),  or 

oo 

m(v0  -  vf)  -  AE  J  tdt  (Eq.  h) 

0 

in  which  v0  and  v/  are  the  (algebraic)  velocities  and  can  have  positive  or  negative 
values  depending  on  whether  the  motion  is  in  the  direction  of  the  intial  impact  or  its 
reverse,  respectively. 

In  general  the  actual  measurement  of  the  pulse  consists  of  a  strain  gage  voltage 
reading  as  a  function  of  time.  This  strain  is  presumably  proportional  to  the  strain 
averaged  over  the  cross  section.  However,  for  reasons  already  discussed,  a  strain  at 
the  gage  itself  may  depart  from  the  average  strain.  Assuming  the  average  planar  strain 
e  is  simply  proportional  to  the  observed  strain  cg  at  the  gage,  i.e.,  «  =  p  eg.  we  can 
write  for  the  correction  factor  f  ""mg  Eq.  (h) 

p  -  m(v0  ~  vf)/(AE  f  Cgdt)  (Eq.  i) 

O 

The  factor  p  ensures  that  conservation  of  momentum  is  preserved.  Wherever  epsilon 
appears  in  the  preceding  equations,  this  should  be  replaced  by  pcg.  Thus,  in  principle, 
the  striker  velocity  z,  penetration  velocity  x,  penetration  depth,  resisting  force,  and 
other  related  quantities  can  be  computed  from  the  impact  pulse. 

Energy  is  distributed  within  the  specimen  and  the  bar  in  several  ways.  Some  of  the 
energy  is  localized  near  the  site  of  the  impact,  e.g.,  that  expended  in  the  non- 
recoverable  (anelastic,  plastic,  viscous)  deformation  process.  There  is  also  a  local  elas¬ 
tic  energy  stored  in  the  impact  region  that  subsequently  is  transferred  back  to  the 
striker  when  it  rebounds.  Finally,  there  is  both  a  distributed  kinetic  and  potential  en¬ 
ergy  associated  with  the  pulse  as  it  moves  along  the  impact  bar.  This  later  energy  can 
be  computed  as  it  develops  with  time.  The  elastic  potential  energy  is  given  by 

tfi-J  (strain  energy/unit  volume)  d  (volume) 

-  1  AEcj  *2dt  (Eq.  j) 

From  the  principle  of  equipartition  of  energy,  this  also  equals  the  kinetic  energy  of  the 
pulse  so  that  the  aggregate  distributed  energy  associated  with  the  pulse  is 

U-  AEcf  t2dt  (Eq.  j') 

0 

This  energy  is  trapped  in  the  bar,  and  excites  the  bar  vibrationally,  ultimately  ending 
up  a  heat 

The  kinetic  energy  of  the  striker  can  be  determined  as  a  function  of  time.  How¬ 
ever,  we  have  found  no  simple  way  as  yet  of  computing  the  amount  of  recoverable 
energy  as  a  function  of  time  during  the  course  of  the  collision.  The  final  irreversible 
work  performed  on  the  specimen,  of  course,  is  simply  the  initial  kinetic  energy  of  the 
striker,  less  its  recoil  kinetic  energy  and  the  pulse  energy  (/. 

Determination  of  the  various  penetration,  velocity,  energy,  and  correction  quanti¬ 
ties  requires  evaluation  of  the  integrals  f  edt,  J  J  edtdt ;  /  t2dt  Once  the  strain  is 
available  in  digitized  form,  the  evaluations  are  straightforward  using  standard  com¬ 
puter  techniques.  The  uncertainties  arise  primarily  from  the  impact  pulse  itself.  In 


principle,  the  pulse  should  rise  from  a  base  line,  representing  the  zero  strain  condition 
and  should  then  return  to  the  same  base  when  the  striker,  upon  rebound,  breaks  con¬ 
tact  with  the  target.  However,  this  idealized  situation  is  usually  not  observed.  Pro¬ 
cedures  for  the  non-ideal  real  pulses  were  discussed  in  5.4. 


■j 

4 


Fig.  4A.  Schematic  representation  of  Impact  bar. 
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Fig.  4.3E.  False  due  to  impact  oi  4.5  mm  steel  ball  with  PC  target  at 
124.5  m/s  as  detected  at  15  cm  gage,  Run  1. 


Fig.  4.3F.  Pulse  due  to  impact  of  4.5  mm  steel  ball  with  PMMA  tar¬ 
get  at  219.9  m/s  as  detected  at  15  cm  gage,  Run  1. 
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Fig.  5.1D.  Plot  showing  lack  of  linear  behavior  of  2/3  power  of  strain  vs  time  In  a 
signal  due  to  impact  of  a  4.5  mm  diameter,  42  g  striker  against  a  PC 
specimen  on  a  PMMA  bar  at  2.1  m/s  P08270C1  and  P08270C2  are 
replicate  runs. 
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Fig.  5. IE.  Plot  showing  lack  of  linear  behavior  of  2/3  power  of  strain  vs  time  in  a 
signal  due  to  impact  of  a  4.S  mm  diameter,  42  g  striker  against  a 
PMMA  specimen  on  a  PMMA  bar  of  2.1  m/s.  P08270P1  and  P08270P2 
are  replicate  runs. 
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Fig,  5.1H.  Plot  showing  linear  behavior  of  1/4  power  of  strain  vs  time  due  to  im¬ 
pact  of  a  4.5  mm,  .37  g  steel  ball  against  a  PC  specimen  on  a  PMMA 
bar  as  detected  at  15  cm  gage.  Impact  velocity  in  12190C1A  was 
124.5  m/s.  Impact  velocity  in  12190C4A  was  219.9  m/s. 
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Fig.  5.1J.  Plot  showing  linear  behavior  of  1/6  power  of  strain  vs  time  due  to  im¬ 
pact  of  a  4.5  mm,  .37  g  steel  ball  against  a  PC  specimen  on  PMMA  bar, 
as  detected  at  30  cm  gage.  Impact  velocity  in  12190C1Z  was  124.5  m/s. 
Impact  velocity  in  12190C4Z  was  219.9  m/s. 
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Fig.  5. IK.  Plot  showing  linear  behavior  of  1/6  power  of  strain  vs  time  due  to  im¬ 
pact  of  a  4.S  mm  .37  g  steel  ball  against  a  PMMA  specimen  on  a 
PMMA  bar,  as  detected  at  30  cm  gage.  Impact  velocity  in  12180P1Z 
was  126.2  m/s.  Impact  velocity  in  12190C4A  was  218.7  m/s. 
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Fig.  S.2A.  Schematic  representation  of  pulse  details  as  observed  at  15  and  30  cut 
strain  gage  positions. 


■ 


Fig.  5.3A.  Difference  between  normalized  actual  strain  and  strain  from  empirical 
equation  5.2C  for  representation  case  of  220  m/s  impact  with  a  PMMA 
specimen. 
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Fig.  S.SA.  Computed  translation  of  pulse  from  30  cm  to  15  cm  gage  position,  for 
case  of  221.5  m/s  impact  with  PC  specimen  pulse  after  translation  is 
shown  compared  pulse  observed  at  15  cm  gage  position. 
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Fia.  6. 1C.  (a)  Corrected  strain  vs  time  for  Indicated  impact, 

(b)  penetration  depth  vs  time  derived  from  above  pulse, 
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Fig.  6. ID.  (a)  Force  vs  depth  as  derived  from  preceding  corrected  pulse, 

(b)  local  energy  absorbed  by  impact  vs  depth  as  derived  from  preceding 
corrected  pulse. 
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Fig.  6. IE.  (a)  Corrected  strain  vs  time  for  indicated  impact, 

(b)  penetration  depth  vs  time  derived  from  above  pulse. 
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Fig.  6.1G.  Force  vs  penetration  depth  for  PMMA  as  a  function  of  impact  ve¬ 
locity  and  type  of  measurement. 
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6.2B.  (a)  Force  vs  depth  as  derived  from  preceding  corrected  pulse, 

(b)  local  energy  absorbed  by  impact  vs  depth  as  derived  from  preceding 
corrected  pulse. 
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Fig.  6.2D.  (a)  Force  vs  depth  as  derived  from  preceding  corrected  pulse, 

(b)  local  energy  absorbed  by  impact  vs  depth  as  derived  from  preceding 
corrected  pulse. 
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Fig.  6.2G.  Force  vs  penetration  depth  for  PC  as  a  function  of  impact  velocity 
and  type  of  measurement. 
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Fig.  7.4A.  (a)  Comparison  of  impact  pulse  synthesized  from  model  with  experi¬ 
mental  pulse  for  220  m/s  impact  against  PC, 

(b)  penetration  depth  vs  time  from  above  pulses. 
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(b) 

Fig.  7.4C.  («)  Comparison  of  Impact  pulse  synthesized  from  model  with  experl 
mental  pulse  for  220  m/s  Impact  against  PMMA, 

(b)  penetration  depth  vs  time  from  above  pulses. 
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Fig.  7.4D.  (a)  Force  vs  depth  derived  from  preceding  pulses, 

(b)  local  energy  absorbed  by  impact  vs  depth  as  derived  from  preceding 
pulses. 
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